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ABSTRACT 


The study undertaken here deals with the deformability of 
rock masses. To attain reliable conclusions on this subject surface 
and underground structures like concrete dams and underground power- 
houses will be analysed. 

At the investigation stage of the construction of large pro- 
jects, like those mentioned above, in situ tests are performed on 
the rock to evaluate the deformability characteristics of a mass. 

The structures we refer to, load a very large volume of rock. Thus, 
large in situ static tests are the ones which will give the most 
representative modulus of deformation for a rock mass. But, how 
representative are the values obtained from an in situ static test 
vis-a-vis the value needed to account for the actual behaviour of 

a structure? An answer to this question is pursued here by under- 
taking the study of case histories. 

For the concrete dams analysed, it is found that the overal] 
modulus of deformation of the rock mass, that accounts for the measured 
structure behaviour, is 1.3 to 5.3 times greater than the static in 
situ test value. For the powerhouse analysed, these ratios vary from 
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CHAPTER I 
INTRODUCTION 


1.1 Rock Mechanics in Engineering 

In the past three decedes, with the coming of larae structures 
like minina excavations, dams, underground aalleries etc., a need 
was felt to study in a more ricorous manner the properties of rock 
as a material and the properties of rock masses as well. The early 
developments of this newly born science, Rock Mechanics, were of 
course areatly influenced by Soil Mechanics which had already pro- 
aressed substantially. Some unfortunate catastrophes made engineers 
look deeper, in particular, into the properties of rock masses. 

The importance of the projects involving the construction of some of 
the structures mentioned above led to the development of better 
instrumentation to characterize the rock mass. 

Rock is the subject of studies from many different areas of 
science. This has brought enaineers, ageo ogists, mineralogists and 
qeophysicists toaether. Today, it is recoanized that to attempt to 
understand the peculiar behaviour of rock masses upon loadina or 
unloadina, the civil enaineer has to be in close contact with the 
structural aeolooist. 

The first two conaresses of the INTERNATIONAL SOCIETY OF 
ROCK MECHANICS have contributed areatly in publicizina answers 
to different facets of the subject and they have also proved the 


need for pursuing the development of theories and tools that will 
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enable us to predict and understand the response of engineering 


Structures involvina the behaviour of rock masses. 


1.2 Behaviour of Engineering Structures 

The study of che behaviour of enaineerina structures can 
be pursued only if a good observation proaram has been undertaken. 
In turn, the quality of that observation program, on a long term 
basis, iS wholly justifiable by the followina three aspects: 

i) Instrumentation will obviously strengthen the safe- 
aquards for public safety, 

ii) Instrumentation will lead to a better understandina 
and improved desicn. A good observation program will often pro- 
vide answers to possible anomalies in the behaviour of a structure. 
The experience oained by enaineers over the years will undoubtedly 
enable them to avoid inefficient construction techniques and tend 
to a better optimization of the cost of the structure, 

iii) The scientific aspect, which of course is indispensable 
if a better understanding of the behaviour of the rock mass and 
consequently, the behaviour of the structure, has to be attained. 

In the following study, we are particularly interested by 
the scientific aspect which is not dissociated from the other two. 
It provides ways of judging the accuracy of the desian hypotheses 
by interpretina the actual behaviour of the structure. 

On the other hand, if the justification of a aood observation 
proaram for lona term conditions is made, it is certainly equally 
important, if not more, to follow very closely the behaviour of a 
structure, durina and shortly after construction. It is in fact 


more important than the lona term observation since at that time, 
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no experience has been aained in evaluating the limits of confidence 
with respect to the initial assumptions made. 

The easiest way of obtainina an overall picture of the 
behaviour of a structure is by measurina displacements. This 
parameter is directly related to the value of the modulus of de- 
formation determined at the design stage. It is the purpose of the 
present study to evaluate the accuracy of the modulus of deformation 
of the rock mass, determined from in situ static tests, in accountina 
for the displacements measured on the structure. Therefore, dis- 
placements will be our standpoint to judge if, whether or not. 
moduli of deformation determined by in situ static tests are 
capable of predicting the correct displacements and, if not, suaqaest 
possible reasons for this discrepancy. 

In this work two particular types of enaineering structures 
will be analyzed i.e. straiaht concrete gravity dams and under- 
around powerhouses. It is obvious that a study of this type has to 
be statistical in that analyses of many structures of a same type 
should be performed if ceneral conclusions are aimed at. By doina 
so, we actually include a variation of the variables affecting the 
behaviour of dams or powerhouses and try to find the common de- 


nominator to the behaviour of those structures. 


1.3 Complexity of the Problem 


The complexity which arises every time that an analysis 
of a structure involvina rock has to be performed is a function of 
two major factors: 

i) the nature of the rock as a material and as a mass, 


ii) our difficulty in taking into account the influence 
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of faults, joints, schistosity etc.. on the behaviour of structures. 

Very often rock shows anisotropy either on a microscopic 
or macroscopic scale. Rock is non-homoaeneous. Rock behaviour is 
time-dependent. 

Up to 1960, no numerical method dared taking into account 
all or even some of the factors mentioned above. The mathematical 
solution would have been too complicated and difficult to apply 
by the average engineer. 

The finite element method entered general use in 1960 
after the work of Clough (1960). This very powerful method has 
increased in popularity ever since because of its versatility and 
capability of handling such problems as non-linear behaviour, 
anisotropy, non-homogeneity, time-dependency and many more. It 
is this method which will be used in the analysis of the dams and 


powerhouses considered here. 


1.4 The Finite Element Method 

It is not the purpose of this section to derive the method. 
Ample literature exists on this subject and may be referred to if 
information is needed. Some references are aiven in the bibliography. 

The structure which is to be analyzed is divided into a 
certain number of elements, interconnected at their nodes. Pro- 
perties are then assigned to these elements. Equations are estab- 
lished for every node in the mesh and each equation includes the 
influence of all the elements on a particular displacement at one 
node. Consequently, a set of simultaneous equations is formed, in 
terms of the nodal displacements, and represents the behaviour of 


the whole structure. The solution of these equations requires the 
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use of a computer. Once the displacements are determined for every 
node in the mesh, it is then possible to evaluate how much distortion 
or strain a particular element has underaone. This is done for all 
the elements. From there, knowing the stress-strain behaviour of 

the material or having made such an assumption, the conversion from 
strains to stresses is made for every element. We therefore obtain 

a complete picture, in terms of displacements and stresses, of the 
structure. 

A finite element analysis can be performed either in two 
or three dimensions. We have said earlier that only straiaht con- 
crete aravity dams and underaround openinas were aoina to be analyzed. 
Since these structures usually have a lenath biager than their other 
dimensions, two-dimensional analyses can be performed on them. By 
doing so, many more analyses can be done which help us in evaluatina 
the influence of various parameters. 

Only linear elastic finite element analyses will be performed 
through this study. Generally, isotropy and homogeneity will be 
assumed. The reasons for working within the framework of linear 
elasticity are the followina: 

i) The conversion of in situ measurements of displacements, 
when doina a static test, to in situ modulus of deformation, is 
commonly done by using the theory of elasticity, 

ii) It is conventional practice in Rock Mechanics to make 
use of linear elasticity in desian, 
iii) A linear elastic analysis is the simplest type of 
analysis one can use and if positive conclusions can be drawn from 


it, it then provides a very easy tool to use for further work, 
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CHAPTER II 


DETERMINATION OF REPRESENTATIVE 
ROCK PROPERTIES 
2.1 Introduction 

It has been recognized, in the field of Rock Mechanics, 
that properties determined on fost samples in the laboratory usually 
overestimate the quality of the rock mass. The reason for this is 
the fact that the discontinuities of the rock mass are not represented 
in the core sample. Therefore, the only logical way of obtaining 
representative rock properties is to go to the field and measure 
them in place. Even then, the problem of finding representative 
properties is not entirely solved. Investigations of the behaviour 
of tne rock mass itself are still needed. 

The deformability of a rock mass is a function of the quality 
of the rock itself and of the characteristics of the discontinuities 
which will be affected by the load transmitted by the structure. In 
the case of a massive rock body, where discontinuities are almost 
absent, the properties of the intact rock will govern the deform- 
ation of the structure. When the opposite happens, for example a 
highly jointed rock mass, the joints will play the dominant role 
in the mode of deformation of the structure. Therefore, since rock 
masses are, in general, fractured determination of joint spacing, 
type of discontinuity, degree of weathering, etc. are of primary 


importance to the engineer who has to decide in which way represent- 
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ative rock deformability characteristics can best be obtained. 

In every project of relative importance, a drilling program 
is included. After the drilling, logging of the cores and subsequent 
interpretation will lead to an overall picture of the quality of the 
rock mass. Deere (1966) has developed a method to describe quanti- 
tatively the quality of a rock mass. This method makes use of all 
the core borings done at the investigation stage. The method is 
called Rock Quality Designation (RQD) and consists of measurina the 
total length of all unweathered pieces of core greater than or equal 
to 10 cm. (4 inches) and dividing the total by the length of the 
particular core run. This index represents an average quality of 
a rock mass and other properties might be related to it. It is with 
this idea in mind that Deere et al. (1969) have investigated the 
possibility of the existence of such relationships. All the dis- 
continuities of a rock mass will affect its deformability. The 
RQD being a measure of the continuity of the body, the above 
investigators have shown that a correlation between the modulus of 
deformation and the RQD existed. Consequently determination of 
the RQD will help estimating the competence of a rock mass. 

Different tests can be performed on the rock to determine 
its deformability. Since static loading represents more closely the 
way rock masses are loaded by the structures we are studying, we 
will concern ourselves with in situ static tests. It is also legitimate 
to think that they are the ones which will give us the most represent- 
ative modulus of deformation for the rock mass. We are not ex- 
cluding dynamic tests from a testing program. They may very well be 


performed and their results can certainly be very helpful in zoning a 
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rock mass. 

In situ static tests are expensive and time-consuming. For 
these reasons the determination of the site where the test will be 
performed is very important. All of the above tools mentioned, like 
core logging, RQD, dynamic tests will help the engineer in deciding 
which volume of the rock mass is likely to have the most critical 
deformability characteristics with respect to the stability of the 
Seructure: 

In the following, a description of the most important in 
Situ static tests and their evaluation will be given. Also, it is 
to be noted that the same terminology used by Deere et al. (1969) 
will be used all through this work. The modulus of deformation wil] 
be the one represented by Ey on the load-deformation curve shown 
in Fig. 2.1.1. This modulus corresponds to the maximum load and 
total deformation. The modulus of elasticity will be the one shown 


as E,_, on the unloading curve. 


2.2 The Different Tests Used and their Evaluation 

For the past twenty years a great deal of work has been 
devoted to increasing the quality of in situ static tests. In the 
following, we will describe plate jacking tests, pressure chamber 
tests, radial jacking tests, thin flatjack test and borehole deform- 
ation tests. The first two have been performed almost systematically 
ever since the undertaking of large projects. The radial jacking 
test constitutes a second phase in the evolution of deformability 
testing equipment. Finally, the thin flatjack test seems to be the 
latest development with respect to large in situ static test. The 


borehole deformation tests will also be looked at. 
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2.2.1 The Plate Jacking Test 


The plate jacking test is certainly one of the oldest methods 
for determining the deformability of rock. We have to distinguish 
here between two different ways of performing this test. The first one 
uses a rigid plate and the second one a flexible pad to apply the load 
to the rock mass. Kujundzic (1955) calls the former a punch test or 
‘uniform deformation test and Serafim (1964) calls the latter a uniform 
pressure test. The punch test has almost completely disappeared in 
current practice due to the complex state of stress it induces ir, 
the rock mass and to the relatively small volume of rock affected by 
the test. Serafim (1966) says that punch tests seem adequate when 
performed on a non-fractured rock mass or on an extremely fragmented 
or very altered rock mass. Since geologic conditions of most rock 
masses are more likely to be within these limits, the uniform pressure 
test is undoubtedly preferable. Rocha (1955) has inserted between 
the reaction system of the plate jacking test and the mortar pad 
applied on the rock surface, oil-filled metallic cushions. This set- 
up is shown on Fig. 2.2.1. For the rock foundation testing of 
Dworshak dam, Shannon and Wilson (1964) have developed a uniaxial 
apparatus which is worth elaborating on. The set-up is shown on 
Fig. 2.2.2. This apparatus included some very special features. 
Firstly, Freyssinet jacks were used to load the rock mass. At the 
same time that loading was performed, uniform pressure on the mortar 
pad was insured due to the flexibility of the jack itself. Secondly, 
since information is needed about the creep behaviour of the rock, a 
pressure recorder-controller device was fixed to the flatjack. As 


the rock was deforming, constant pressure in the jack was maintained 


of 
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by the use of a hydraulic pump. Thirdly, the way deformations were 
measured is of particular importance. As shown in Fig. 2.2.2 AX size 
boreholes were drilled at both extremities of the apparatus and Carlson 
jointmeters were installed in them. The jointmeter usually extended 
to a depth of 5.5 m. (18 feet) from the face of the gallery and con- 
sequently allowed for the measurement of an average deformation for 
a large mass of rock. The usual diametral extensometer was also 
used. 

The U.S. Bureau of Reclamation had at the same time developed 
a very similar apparatus as the Shannon and Wilson one. Wallace et 
al. (1969) give a good description of it. The main difference is in 
the way deformations were measured. The principle is the same but 
more information is obtained by the U.S. Bureau test. Rock deform- 
ations were measured by a diametral extensometer and two seven 
positions retrievable extensometers (REX-7P). The diametral exten- 
someter, as well as the REX-7P were installed directly underneath 
the loaded surface. For this arrangement, small holes were provided 
in the middle of each flatjack. The retrievable extensometers were 
fixed in a 6 meters long NX hole and deformations were measured 
between every one of them and the collar of the hole. The measuring 
device consisted of seven linear variable differential transformers 
(LVDT). This way of measuring deformations has proved to be extremely 
useful. 

It is known that whenever a gallery is excavated, the initial 
state of stress in the rock mass is disturbed. This produces a so- 
called “destressed zone" (Deere et al., 1969) around the opening which 


influences greatly the deformability measured by an in situ test. 








gxt2 KA $.8.8° -pET af inwotie 2A al oe 
noz fed bos 2uteysqqs att to agtanaenap alia te belfivb sew esfortared 
babiotxe yf favey vet amIhiot srt) wo 8a aent ane eratomintos, 
-noo bee yvallép Sit Yo e067 ong mot (Yaar Bi), a2 to digqeb 6 oF 
+10} nottamvoteb sperevs ns 76 dnomowesom odd “ar bawol ls xltnoupse 
o2ls 2aw 19dsmoO2nstx9 Tergemste fau2d nT door to 226m apyet ae ) 
-b92u — 

bsqofsvab. srt omse ot ts bad notes foah Fo yssqyud 1.2.0 onT : 
ta sosllsW .sno noeliW base nonnsn2 sit 26 euservéqas yelimre yvev 6 
at 2t sa9netsttib ntem oAT tr to nottqtyaesb boop s svkp (ROOT) fs 
sud sms2 sdt et ofqtoniyg anil _bswebom 919w anol tsmvoteb Yow ans 7 : 
-myoteb As0h .teat usowé .2,l) sat yd bontstdo 2f nottsmiotnr Sion 
havez ows bis yatemoe2natxs fevtoms th 6 YU bo wW269M s1ow anette 
-n9dx5 [sudematb oAT .(9%-X39) erat emoanetxa sTdeverytey 2notitzog 
atesnvsbay yfrosytb bel stent ate qy-N3A ony es [isw 26 ,veramoe 

bebivorg a1aw eofor [fame eT RSMSP 616 2ind 403 .poRtwe babsol ons a 

evsw 2vetemoenstxs sidevaiatsy adT | aebielt ines to afb tm and ore” @ 

bsrwe5sm svaw enol temrotsb bas slot XK paof evotam a » at bextt 
pniwwesem sdT .9for sit to Yel Too ont brs wart to sno \vove nsewisd : 


evonotensys [sidnsastt tb ofdsivay agent | nevse To beselenoo solveb on 
ene ad ot bovorg asi 2noitsmoteb pniiwesam to Di 2taT (roy) : 





14 


Benson et al. (1969) quote for one test in particular that 33% of 
the total deformation between 0 and 4.6 meters (15 feet) from the 
rock surface had occurred in the first meter. Therefore, it is 
advisable to measure deformations in depth from the rock surface in 
a test adit. This will obviously lead to a more representative value 
of the modulus of deformation of the rock mass. 

In brief, it is suggested that the best quality in situ 
plate jacking test that can be performed now on a rock mass is one 
of the same type as the U.S. Bureau of eration test. The 
aluminum tubing makes the apparatus light enough for easy trans- 
portation; the Freyssinet flatjacks arrangement fulfills a triple 
purpose and finally the deformation measuring device provides a 


better picture of the rock mass behaviour. 


2.2.2 The Pressure Chamber Test 

The pressure chamber test is one of the largest in situ tests 
performed on a rock mass. For this reason, this test is used much 
less frequently than the other in situ static tests. However, the 
importance of a project together with the rock conditions at a site 
can easily justify its use. Loading a very big volume of rock, this 
test will be of particular interest when a large joint spacing 
characterizes the rock mass to be tested. Figure 2.2.3 gives the 
usual set-up for the pressure chamber test. 

Application of the pressure on the rock mass is achieved by 
filling the chamber with water. From this comes one of the most 
delicate problems of the test, to assure the water tightness of 
the chamber. In some of the earlier pressure chamber tests performed 


(Kujundzic, 1955), water pressure was directly applied to the rock 
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surface. The amount of water collected after the test was then used 
to estimate the permeability of the rock. Rocha (1955) has used 
bituminous material as an impermeable liner. Italian engineers seem 
to be the ones who have used this test the most extensively. Figure 
2.2.3 was based on the usual set-up used in Italy. As seen on this 
figure, water leakage was prevented by using a rubber sealing bag. 

A concrete lining with transverse and longitudinal joints was also 
used between the rock surface and the sealing bag. 

In the pressure chamber test, deformations are often 
measured at the middle section of the chamber. This has for effect 
to minimize the influence of the extremeties of the loaded area. 

The deformation measuring system usually consists of extensometers 
installed on four different diameters, one being vertical, another 
horizontal and the others at 45° to the first two. Usage of buried 
deformation gages, for the pressure chamber test, hasnot been found 
in the literature but it is expected that they have been used in 
the latest tests. The advantages of using these gages have been 
-discussed previously. 

To summarize, if provision is made for measurement of 
deformations at depth, from the rock surface, the pressure chamber 
test is likely to give the most representative rock modulus of 


deformation. 


2.2.3 the Radial Jacking Test 


The radial jacking test is a combination of a plate jacking 
test with a pressure chamber test. On one side, it uses the relatively 
simple mounting procedure of the plate jacking test and on the other 


benefits of the large volume of rock loaded by the pressure chamber. 
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The first radial jacking apparatus was developed by the 
Hydrotechnical Institute of Belgrad (see Kujundzic, 1955). Although 
almost twenty years old, the test had much similarity with the one 
performed today. Figure 2.2.4 shows schematically the different 
characteristics of the test set-up. As seen on the figure, the 
gallery had a diameter of 2.20 meters. The load was applied on a 
concrete lining by sixteen flatjacks having a length of 1.85 m. 
each. The inside face of the jacks was fixed to cylindrical wooden 
reaction system. Transmission of pressure to the flatjacks was 
obtained by manually pumping water into them. The deformations were 
measured at five different sections, three of which were directly 
under the influence of the flatjacks. At each section, the deform- 
ation measuring device consisted of four extensometers making 45° 
angles between them. Finally the interpretation of the deform- 
ation results accounted for the fact that the test had a small 
LENGTH/DIAMETER ratio. 

The radial jacking test has also been used extensively, 
for the design of pressure tunnels, by Tiroler Wasserkraftwerke, 

A. G. (TIWAG) in Austria. The apparatus they have used was basically 
equivalent to the Yugoslavian one but a few modifications were needed. 
Lauffer and Seeber (1961) give a short description of the Austrian 
apparatus. Since high internal pressure can be applied on the rock 
surface in pressure tunnels, the main changes consisted of re- 
designing the reaction system of the apparatus. Strong hardwood 
planks and ring supports of high-tempered steel were used. Those 
modifications then allowed for the application of stresses in the 


order of 60 Kg/cm°. The other features of the test with respect 
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to the ones described by Kujundzic (1955) seem to be very similar. 

The U.S. Bureau of Reclamation has also developed its own 
radial jacking test. Wallace et al. (1969) describe the Bureau 
test. Again, it resembles the first two tests mentioned above with 
one very important difference: deformations are measured in depth 
from the rock surface to eliminate the effect of the destressed 
zone surrounding the excavated gallery. Measurement of diametral 
deformations of the gallery is also provided. This way, depending 
on the problem which has to be solved, i.e. the design of a pressure 
tunnel or a surface structure, the corresponding deformations are 
used to evaluate the modulus of deformation of the rock. 

The radial jacking test, as it is performed today, is 
certainly the best optimization one could have attained from the 
plate jacking test and the pressure chamber test. As we said 
earlier, the assembly is relatively simple; the test loads a large 
volume of rock; the test area is accessible to the operators; al] 
of the equipment is recoverable. On the other hand, the cost of 


such a test is high. 


2.2.4 The Thin Flatjack Test 
The thin flatjack test has been developed mainly by the 


Laboratorio Nacional de Engenharia Civil (LNEC). The procedure used 
in performing the test is as follows (Rocha et al. 1970). 
i) a diamond disc saw of diameter d is used to open one or a 
group of slots, contiguous and in the same plane perpendicular to 
the rock surface, having a depth h and a thickness t, 
ii) thin flatjacks of height (h - a) are then introduced in 


each slot, 
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iii) a pressure p is applied to the flatjacks, 
iv) deformation of the slot is measured by four deformeters 
inserted in each flatjack, 

v) reduction of deformations to modulus value is made 
through diagrams. 

A representation of three thin flatjacks inserted in the 
VOCK iss aivenein) Fide 2.c-5. 

Like the other in situ static tests, the thin flatjack 
test will vary with the nature of the rock and the spacing of the 
joints. Rocha (1969) has shown the relative decrease of the modulus 
of deformation by extending the length of the slot. 

In general, the thin flatjack test has proved to be of a 
good quality to determine the modulus of deformation of a rock mass. 
The apparatus can be taken in a gallery and can measure represent- 
ative properties providing the slot is cut to a depth beyond the des- 
tressed zone. Preparing and performing the test is done relatively 
quickly and at a much lower cost than the first three tests men- 
tioned previously. 

It is felt that the important weakness of this method is in 
the interpretation of the deformation results. Rocha et al. (1970) 
have assumed that the rock mass was cracked in the plane of the slot 
thus allowing them to represent the test by a quarter-space under 
pressure p. They have developed a special technique to interpret 
the thin flatjack results and suggested that care should be taken 
when performing the test in a region of high compressive stresses 
perpendicular to the slot. It would then be advisable to enlarge 


the size of the slot all-around the proposed loading area. 
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redeks Borehole Deformation Tests 

When the size of a project will be such that the use of a 
large in situ static test will not be justified, borehole deformation 
tests may be used for determination of the modulus of deformation. 
Table 2.2.1 enumerates some of the borehole tests used and gives 
some characteristics of each. 

As can be seen, the borehole deformation tests are divided 
into two major groups depending upon whether transmission of load 
to the rock is done through a rigid or flexible component. The main 
differences between the two groups are the following: 

i) Usually a higher stress level will be attainable when 
using a rigid apparatus, 

ii) the non-uniformity of the state of stress induced by 
the rigid test into the rock mass will be more pronounced than for 
the flexible test, 

iii) although possible in the case of a flexible test, 
cracking of the rock will be more probable when using a rigid 
apparatus. 

iv) for equal size tests and for an equal average applied 
stress, a smaller volume of rock will be loaded by the rigid apparatus 
than by the flexible one. 

The first two items, will be partly responsible for crack- 
ing of the rock. Moreover, the theory of elasticity says that a 
radially applied pressure in a hole will induce a tensile tangential 
state of Stress jaround, 16. 

Cracking around the borehole for a flexible test also has 


to be included in the interpretation of the results. Rocha et al. 
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(1966) have discussed of cracking around the borehole and have 
included its effect in their evaluation of the modulus of deform- 
ation. 

With respect to the volume of rock tested, we have said 
that the flexible borehole test would load a larger volume than the 
rigid one. Even then, interpretation and experience will tell if 
whether or not the rock conditions involved in the flexible test are 
representative of the whole rock mass. Extreme care has to be taken 
in that direction since totally erroneous values could be used for 
design purposes. 

The method of drilling the borehole, especially in soft 
rocks, may influence the result. Meigh and Greenland (1965) give 
some values of the modulus E measured with the Menard pressure- 
meter in a sandstone deposit. For the drilling of one Porenote! 

a non-coring water-flush technique was used. For two other bore- 
holes, a drag bit with air-flush and a thin-walled diamond coring 
bit with water-flush were used. The results presented show that 
the moduli measured in the boreholes drilled with the non-coring 
water-flush technique are consistently lower than when using 
either one of the other methods mentioned. That is to say, the 
substantial amount of water used during drilling, in the first 
case, has probably resulted in some softening of the material 


around the hole. 
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Cus) “DETaTl| Se Oreiint@s 1 tum reste Data 
and Comparison of Results 


Up to 1966, data on static in situ tests have been collected 
and presented in a report by Deere et al. (1969). This has been ex- 
tended here to 1973. The data are presented in Fig. 2.3.1. Comparison 
of the modulus of deformation and the modulus of elasticity is made in 
that figure. As can be seen the newly collected data agree relatively 
well with the line suggested by Deere et al. 

In the same report mentioned above, the investigators have 
compared some in situ moduli to an overall modulus of deformation com- 
puted from displacement observations of structures. In the period 
going from 1967 to 1973, only two cases were found where the same 
comparison could be made. Figure 2.3.2 shows the points plotted by 
Deere et al. (1969) together with the earlier ones. Even if only 
a few points are plotted on that graph, a very interesting feature 
emerges. Most of the points representing static in situ tests, 
except two, are lower than the 1:1 line. This means that if de- 
signers have used those moduli values, determined by the static in 
situ tests, a conservative design has resulted. No comment has been 
given by Deere et al. with respect to the highest point which cor- 
responds to a plate jacking test. On the other hand, Kleiner and 
Acker (1971) describe the different characteristics of the rock 
mass at Mossyrock dam site in Washington state, U. S. A. The point 
on the 1:1 line corresponds to a plate jacking test performed on a 
massive andesite. As suggested by the authors, the agreement of the 
modulus determined by the large scale in situ static test with the 
One computed from the displacement measurements, is probably due to 


the natural massiveness of the rock body of Unit B. 








-Ka nead esd 2taT _(eaer) “fs t9 seat xd ne 6 nt bednezsyq bas 
noetyaqmod .[.€.S .pra nF besyns2ay¢ S16: 546 sat BNET os ayad bebnes 





at abem 2¢ yitortesfa to eufubom srl bn6 nofsenrrotab. to zulubom ont to 
vlavitslay osips stsb botoslfoo \ Iwan ant nase ad, ie) 2A ..ouuntt ted? 


.fs 5 59980 yd bet zene ant ont astw Tisw 


aver aiodspit2avat edt .svods benoftnsm droge ompe sad ni © 

-moo notdanmioteb to euTubom [lsveve ns od rfubom utte nt smoz bervsqmoo 
doiysg oft nl .esiwiouryse to enotisvysedo tnsmezsigerh mort batug 
ams2 sit siedw bavot ovsw 29269 owt yIno .ENeT ot Yael mort entop 


vd bettolq edntog ant eworle $.£.8 sword  .sbem od bluos neeiveqmoa 


vino ti mova .2eno refiyss ons dttw vedtepos (C021) is Je sised 
siudes? onfdesiaint yay s .dqexp Jedt no beddolq evs etntog wat 6 
eet2at utt2 nf otist2 patingesiqey eintog alt to yeOM .2spem— 

~sb ti dedt ansom 2tdT ontl T:f art ned vowol ays end JqaoKe 

nt sftes2 ont yd bsntmryadsb .2sutsy flubom sent bezy svar 2xafipte 
nesd 261 Jnanmoo ol ss vest 26f éeienb ayitavisenoa. 6 e2teet, udke 
~Y02 notriw tatog +eargte ont ot $osq251 fitiw .Is +9 seed yd gavin. 
brie: senrety brent nando, ost nO .test onttost sisiq 6 od ebnagest 
daoy itt Yo eotszitatoersdo Jnovettib ont odtyogsb (TVEL) a eaaA 
‘tntog siT «A .2 .U -atete nosontreslt ni atte me Joowz2oM 46, 226m 
sab betes enisost apes io muchetitie 09 ne et -“, 80" 


on a 
: 4 A a on is ae. wu - 
mee Bs oe om =a Sa) Fa ? |e v- rT 
“5 a out ; 2 ist - , 
: : " ¢ if 
“ A ; . 





Modulus of Deformotion (Field, Static), Eg» psix 10° 





fol] 
0.1 









































0.2 ; 06 08 1.0 2 4 6 8 10 
Modulus of Elasticity, (Field, Stotic), E,, psi x 10° 


PLATE JACKING TEST OQ 
RADIAL JACKING TEST @ 


PLATE JACKING TEST ® 

RADIAL JACKING TEST © 

BOREHOLE DEFORMATION TEST 4 
(UP TO 1969) (1969 to 1972) 


FIG. 2.3.1 COMPARISON OF THE MODULUS OF DEFORMATION AND 
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The points representing the moduli at Tehachapi Tunnel, in 
California, Were compared to a modulus obtained by the tunnel relax- 
ation method. This method is described by Kruse (1969). 

On very few cases a comparison has been made between the in 
Situ static modulus and the one to account for the measured deform- 
ations. For the cases studied previously we see that the behaviour 
of the structures involved correspond to an overall modulus higher 
than the one measured by static in situ tests. This fact supports 


the concept of the study undertaken herein. 


2.4 Determination of Pre-Existing Stresses in a Rock Mass 


Knowledge of the natural state of stress in a rock mass is 
absolutely necessary when dealing with underground openings. Further- 
more, if a better understanding of the behaviour of the earth's crust 
is of interest, determination of the natural state of stress in rock 
masses is mandatory. 

In the past two decades extensive work in developing equip- 
ment capable of measuring the natural state of stress has been carried 
out. Obert (1966) has summarized the principal methods used. New 
developments were presented at the International Symposium on the 
Determination of Stresses in Rock Masses, held at Lisbon in 1969. 

No significant improvement has taken place since. Then, the following 
is based mainly on the above two references. 

In the determination of the natural state of stress in a rock 
mass one has to distinguish between two major groups of methods: 
the direct method, i.e. a method in which a stress is measured and 
the indirect method which makes use of the elastic properties of the 


material to obtain the state of stress. In the first category we have 
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rigid inclusion stress gauges and flatjacks. The second group con- 
sists mainly of displacement measuring devices and strain measuring 
devices. It is not the purpose of this section to describe all of 
these instruments. The two references mentioned above should be 
referred to if detailed information is needed. However, the follow- 
ing will give a short description of one instrument belonging to 


each group. Also discussion of the latest developments is presented. 


2.4.1 The CSIR ‘Doorstopper" 
In 1965 the South African Council for Scientific and Industrial 


Research undertook the development of a rock stress measuring device. 
Commercially available since 1968, the CSIR "Doorstopper" has already 
been used all over the world. 

The “doorstopper" has been designed to determine the absolute 
stress in rock, using an overcoring stress relieving technique. The 
method is therefore an indirect method. Performance of a "door- 
stopper" rock stress measurement is as follows: 

i) a BX borehole is drilled to the required depth and its end 
is flattened and polished with diamond tools, 

ii) with an installing tool the "doorstopper" strain cell is 
bonded on to the end of the borehole and initial strain readings are 
recorded, 

iii) the borehole is then extended with a BX diamond coring 
crown resulting in a stress relieving of the core, 

iv) after overcoring, final strain readings are taken and 
removal of the core and "doorstopper" is done, 

v) knowledge of the elastic constants of the rock, when over- 


cored (see Benson et al., 1970) will allow for the determination of 
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the stresses present on the flat end of the borehole before overcoring, 
vi) a relation exists between the above determined stresses and 

the stresses of the surrounding rock. Use of this relation permits 

the calculation of the surrounding stresses. Different relations 

have been summarized by Leeman (1969). 

Rock stress determination with the "doorstopper" can be done 
in one or three boreholes. If the measurement is performed in one 
borehole, it is implicitly assumed that the borehole axis is parallel 
to a principal stress. Consequently, the flat end of the borehole 
is a principal plane and the problem is reduced to a two-dimensional 
one. In the case where the validity of such an assumption can be 
questioned, three different boreholes need to be drilled in order to 
determine the complete state of stress at a point. 

In item 5 above, we see that the elastic constants of the rock 
are needed. A "hydraulic assimulator" (see Leeman, 1969) has been 
developed with which it is possible to calculate the state of stress 
acting around the flat end of the borehole before being overcored. 
The important feature of using this assimulator is that the elastic 
constants of the rock are not needed. In brief, the "hydraulic 
assimulator" loads the core, with the "doorstopper“ glued on it, in 
such a way as to restore the initial strain readings. The load 
needed to do so is directly related to the principal stresses acting 
around the flat end of the borehole. 

The reliability of rock stress measurements with the “door- 
stopper" seems to have been proved (Heerden, 1968). However the 
stresses measured locally on a small piece of core need not be 


representative of the stresses in the mass at a larger scale. 
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2.4.2 Flatjack 


As mentioned previously, the flatjack technique is referred to 
as a direct method of rock stress measurement. Determination of 
the stresses with the flatjack method is done in a gallery excavated 
in the rock body. As is well known, the excavation of the gallery 
together with blasting will greatly influence the state of stress of 
the rock, in the vicinity of the opening. Therefore, stresses measured 
with the flatjack technique, which is performed closed to the rock 
surface in a gallery, will not represent the state of stress of the 
rock away from the zone of disturbance. As underlined by Rocha 
(1969), measurements of the state of stress close to the rock surface 
by the flatjack method will be applicable to tunnel-driving methods, 
tunnel maintenance and tunnel-lining design and construction. 

In principle, the flatjack method consists of: 

i) fixing references (vibrating wire gages, steel pins, strain 
gages, etc.) on the rock surface of an adit, 

ii) taking initial readings (frequency of the wire, distance 
between pins, initial strains, etc.), 

iii) making a slot in the rock mass, close to the references, 
and inserting a flatjack in the slot, 

iv) pressurizing the flatjack so that the initial readings are 
recovered. 

The pressure necessary to reach "cancellation" is then con- 
sidered as the natural stress of the rock mass, perpendicular to the 
flatjack. The stress found by the flatjack technique is an average 
stress acting over the relatively large area of the flatjack. For 


this reason, the flatjack technique is rather insensitive to local 
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concentration of stresses. 


2.4.3 Recent Developments in In Situ Stress Determination 


In the last ten years, two different organizations have 
developed a technique by which the complete determination of the in 
Situ state of stress could be obtained in a single borehole. LNEC 
worked on a plastic cylinder in which an assembly of electrical strain 
gauges was embedded. The CSIR built a triaxial cell capable of 
measuring strains with three rosettes at a specific relative orien- 
tation. 

Undertaking a measurement with these two methods is basically 
the same. It consists of an overcoring technique. The LNEC plastic 
cylinder can measure initial as well as induced states of stress in 
the rock mass. The CSIR triaxial cell is capable of measuring the 
initial state of stress only. 

Further details are given in Rocha and Silverio (1969) and 


Leeman (1969). 
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CHAPTER III 
ANALYSIS OF CASE HISTORIES 


3.1 Introduetion 
Generally, the selection of case histories in this study was 
based on the following requirements: 
i) results of in situ static deformation tests had to be 
available, 
ii) loading on the dam or in situ rock stresses had to be 
known or estimated, 
iii) the deformation behaviour of the structure considered 
had to be given. 
On this basis three dams and one underground powerhouse 


were selected for analysis. 


3.2 Analysis of Dams 


The three dams that are to be analyzed are Krasnoyarsk dam 


in East Siberia, Alpe Gera dam in Italy and Bhakra dam in India. 


3.2.1 Boundary Conditions 


In the boundary conditions we will consider the mesh dimensions, 
the displacements boundary conditions and the loading boundary con- 
ditions. All of these items are shown in Fig. 3.2.1. 

If H is the height of the dam, the boundaries of the finite 


element mesh were as follows: 
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FIG. 3.2.1 BOUNDARY CONDITIONS USED IN THE ANALYSIS OF A DAM. 
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i) the vertical sides of the mesh were set approximately 
2.5H downstream of the toe and upstream of the heel, 

ii) the bottom boundary of the mesh was about 4.0H below the 
base of the structure. 

On the vertical sides of the mesh horizontal movements were 
prevented and vertical movements were allowed. All along the bottom 
boundary displacements were restrained in all directions. 

The loading boundary conditions consisted of water forces 
applied on the upstream face of the dam, on the surface of the 
reservoir bed and where tailwater was present, on the rock surface 
downstream of the structure. Also the influence of uplift pressures 
has been simulated by applying upward forces to the base of the dam. 
Naturally, this is only approximate. An effective stress analysis 
would take uplift into consideration more accurately. However the 
accuracy of such an analysis would depend on the accuracy and number 
of water pressure measurements in the rock mass. No common rule, 
with respect to the uplift pressure, is set for the three dams that 
will be analysed. The uplift forces applied to the base of every 
structure have been taken from actual measurements. 


All the analyses assumed plane strain conditions. 


3.2.2 Procedure Used in an Analysis 


Depending if the displacements of the structure are given 
with respect to a time zero, corresponding to the beginning of the 
construction of the dam, the procedure used will vary. Such a 
value for displacement will include the effect of gravity and the 
subsequent water loading. Therefore, if for a structure, displace- 


ments (referring to a time zero) are given for one particular time in 
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its history, only one analysis need to be performed. This analysis will 
include gravity loading and the water load acting on the structure 
at that particular time. On the other hand, if, for example, dis- 
placement measurements have been interrupted for a certain period, 
in the history of the structure, a different procedure will apply. 
This procedure consists of performing two analyses, each one represent- 
ing different loading conditions in time on the structure. By doing 
SO, we consider the first analysis as a reference to the second 
analysis and we look at the relative displacement. 

For the analysis of Krasnoyarsk, the first case mentioned above 
applies. Alpe Gera dam and Bhakra dam are analysed for two different 


dates in their history. 
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3.3 Analysis of Krasnoyarsk Dam 


3.3.1 Introduction 

Krasnoyarsk dam is one of a series of six dams planned on the 
Enisei river, in East Siberia, to regulate its flow. The dam also 
serves the purpose of producing electricity. In fact Krasnoyarsk 
station has the highest production in the world. 

The dam is a massive concrete structure with a height of 
124 m. and a base width of about 100 m. It has a triangular profile 
with a vertical upstream face and an inclined downstream face, having 
a slope of 1:0.8 in the spillway section and 1:0.76 in the power 


house section. It is founded on granite. 


3.3.2 Method of Investigation and Values 
Assumed for the Study 


As described by Bochkin et al. (1971), the dam is founded 
on granite that ranges from sparsely jointed to highly jointed. 

The modulus of deformation of the foundation rock has been 
found by subjecting a rigid plate to vertical loading. This test 
has been performed on the different rock conditions mentioned above 


and gave the following results: 


i) for sparsely jointed granite E = 160,000 Kg/em? 
ii) for medium jointed granite E = 90,000 Kg/cm® 
iii) for highly jointed granite E = 45,000 Ka/cm® 


Poisson's ratio, for this granite, was assumed to be equal 
to 0.2 in the finite element analysis. The unit weight of the rock 
was set to zero because only the displacements caused by the change 
of stress in the rock mass is of concern here. 


Two other materials had to be specified in the analysis. 
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Obviously, the concrete of the dam body was one, and the crest 
gate area on the dam was the other. 

The elastic tangent modulus of the concrete was taken equal 
to 360,000 Kg/cm, as determined on 10 x 10 x 40 cm. samples at 
the age of 28 days. As can be seen in Table 3.3.1, the tangent 
modulus for this concrete, more than doubles in a period going from 
three days to 360 days. 

The effect of the stiffening of the concrete in time, on 
the deformations of the rock mass was studied by doing two analyses 
which will be discussed in the following section. 

Poisson's ratio and the unit weight of the concrete were 
assumed to be 0.2 and 2.4 Tm/m? respectively. 

The crest gate area is expected to be more deformable than 
the dam body itself. Since it does not reproduce in the third 
dimension, a plane strain assumption for this section is not exact. 
The plane stress solution can be obtained by changing the elastic 
properties of the concrete, using the theory of elasticity. In 
order to account for this, the modulus and Poisson's ratio for the 
crest gate area have been assumed to be 250,000 Ka/em® and 0.25 
respectively. However, the results proved to be insensitive to 
that change. Finally, since this area of the dam is more porous 
than its body, a weight of 2.0 Tm/m> was assumed in the analysis. 

Different analyses involving different conditions on the 
life of Krasnoyarsk dam are presented in the following section. 

The finite element mesh used through the study is shown in Fig. 
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3.3.3 Results of Predicted and Measured Displacements 


As stated in the previous section, a study of the effect 
of the stiffening of the concrete in time, on the foundation dis- 
placements, was done. Two analyses were performed assuming an E 
value determined after 28 and 360 days. These values are given 
Wtablew3. 371. 

The results of the displacements of the foundation surface 
for the two cases are given in Figures 3.3.2 and 3.3.3. It is 
obvious that, in the range of variation considered here, the fact 
that the concrete stiffens as it ages does not affect the deform- 
ations of the rock surface. In the case of a modulus equal to 
360,000 Ka/em®, the average vertical displacement over the dam- 
foundation contact is 21.1 mm. The second analysis used E = 
470,000 Kg/cm® and an average displacement of 21.0 mm. was cal- 
culated. Since the difference was insignificant, the value of 


360,000 Kg/cm® was assumed in the following analyses. 


i) Gravity Analyses 


Due to the availability of foundation displacements 
under the weight of the dam itself, a gravity analysis was under- 
taken. The modulus of deformation used here was an average of the 
three values mentioned above, i.e. 1.0 x 10° Kg/cm*. The dis- 
placements are plotted and shown in Fig. 3.3.2. The vertical dis- 
placements over the whole dam-foundation contact averaged at 21.1 mm., 
with a maximum value of 23.7 mm. The settlement measured under 
section 22 of the dam was 10.2 mm. 


It can be concluded from these results that, if the modulus 
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of deformation of the rock mass below section 22 is taken as the 
average of the values determined in situ, the displacements pre- 
dicted will be 2.3 times the measured ones. This implies that, 
to account for the measured movements of the foundation, its 
modulus should be equal to 230,000 Kg/cm*. 

Based on a linear elastic analysis, we will always find 
that no matter which modulus is assumed, if we want to match the 
measured settlements, a stiffness of 230,000 Kg/cm® will be arrived 
at., The fact that the end result is always the same is rather 
interesting because the only variable involved here is the decision 
on the design modulus. 

It is very unlikely that the designers of a dam will use 
the maximum value of the modulus of deformation found from in situ 
tests. In the case of Krasnoyarsk dam, even if they had chosen 
the highest value, i.e. 1.6 x 10° Kg/cm*, a difference of about 
fifty percent would have still existed between their choice and 
the stiffness arrived at earlier. This is shown on Fig. 3.3.4 


which summarizes the results of the gravity analyses. 


ii) Gravity and Water Load Analysis 


Filling of the Krasnoyarsk reservoir began in April 1°67. 
The history of the impounding is shown on Fig. 3.3.5. 

The available data on Krasnoyarsk dam for section 22 of 
the spillway is presented for July 1969, right after a steep rise 
in the water level. Readings of the uplift pressures on the dam — 
base have shown that the grout curtain and the drains built in the 


foundation have been very efficient in reducing them. Even if the 
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uplift was very small, i.e. approximately one per cent of the 
total head at that time, uplift forces on the dam base were in- 
cluded in the finite element analysis. 

The properties assigned to the crest gate area and the 
dam body were the same ones as for the gravity analysis. The rock 
deformation modulus was taken as 100,000 Kg/ cm? with a Poisson's 
ratio of 0.2 and no unit weight. 

The loads involved in this first analysis are then, the 
dead weight of the dam, the water pressure acting on the upstream 
face of the structure and on the reservoir floor and lastly, the 
uplift pressures on the dam base. 

The total displacement vectors of the foundation are shown 
schematically on Fig. 3.3.6. The maximum settlement calculated 
was 40.8 mm. and the average over the dam-foundation contact was 
37.1 Ine 

If, for example, we assume that no hydrostatic pressure 
was transmitted to the grout curtain, which is very unlikely and 
also inconsistent with uplift, we would find that, to account for 
a total displacement of 19.4 mm. up to July 1969, the overall 
modulus of the foundation rock should be 210,000 Kg/cm*. The 
value of 19.4 mm. and others are given in Table 3.3.2., from the 
beginning of impounding up to July 1969. The more reasonable 
assumption of full hydrostatic pressure being applied on the grout 
curtain will be looked at later. 

An interesting feature appears if we look at Table 3.3.2 
together with Fig. 3.3.5. Up to July 1969 the two severe loading 


periods are from April to August 1967 and from May to July 1969. 
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The year 1968 could be qualified as one of constant loading. On 
the other hand, a displacement of 3.3 mm. is given for that period. 
This could mean that the rock in the foundation is creeping under 
that constant load. Since the analysis does not include any time- 
dependent calculation of the deformations the 3.3 mm. could be 
subtracted from 19.4 mm., which would bring the overall E to 
approximately 2.5 times the initial assumption. 

The creep effect mentioned above versus the linear elastic 
finite element analysis can be looked at from another point of view. 
We seek a correlation between the numerical analysis and the actual 
deformations measured on the structure. This correlation has to 
embrace all the possible factors contributing to the deformations 
of the dam such as discontinuities, anisotropy, time-dependency. 

We therefore include the possible creep mentioned above in the 
total displacement and keep the ratio of 2.1 between the overal] 


modulus and the one used in the study. 


iii) Gravity and Water Load plus Assumption 


of a Crack in the Foundation Rock 


The final analysis performed on Krasnoyarsk dam is a more 
realistic one. It considers the foundation rock to be cracked to 
a depth equal to the height of the dam and having full hydrostatic 
pressure applied on the grout curtain. The dam is then free to 
move under the varying water levels. It is separated from the 
rock mass upstream of it by an area of zero rigidity. Consideration 


of the hydraulic aspects of this was given by Casagrande (1961). 
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The displacements experienced by the dam upon impounding 
up to a reservoir depth of 105 meters are shown on Fig. 3.3.7. The 
maximum settlement is 31 mm. and the average over the dam width is 
30 mm. If we compare this displacement with the one mentioned 
before i.e. 19.4 mm., we obtain a ratio of 1.6 between the predicted 
and the measured. Therefore, if the conditions involved in the 
analysis are realistic, an overall modulus of 160,000 Kg/em? would 
make the predicted settlements essentially equal to the measured 
ones. 

Finally, the Russian investigators have also found by con- 
Sidering the dam-foundation inclination due only to water loading, 
that, to account for the measured deformations of the rock at section 
22, the overall modulus of deformation had to be equal to Zn xX 10° 
Ka/cm?. This modulus, with respect to our assumptions, gives < 
a ratio of 2.54 which is in reasonable agreement with the values we 
have quoted all through the Krasnoyarsk study. Table 3.3.3 sum- 
marizes the results of the finite element analyses and of the 


actual behaviour. 


3.3.4 Comparison of Actual Behaviour With Predicted 


In the case of the Krasnoyarsk dam, as in the case of any 
other dam, a representation of the actual deformations of the found- 
ation of the structure and of the structure itself will be faced 
with two major alternatives: 

1) either the rock mass on which the structure is sitting 
is initially very sparsely cracked and has little tendency to open 


in tension upon impounding of the reservoir. This is a function in 
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particular, of the initial state of stress, the E pS 


ratio, of the tensile strength of the rock and of the water pressure 


concrete 


distribution in the rock mass. 

2) or the rock mass is initially highly cracked and/or 
does not restrain the development of crack formation at the heel 
of the dam. 

The first case would make the dam body integral with the 
rock foundation, which in turn would contribute, by its stiffness, 
in reducing the deformations of the structure. The second case would 
"disconnect" the dam body from the upstream rock, and make it very 
much easier to deform and rotate. 

If we refer to the results obtained in the analyses of 
Krasnoyarsk dam and assume that, from the beginning of impounding 
there is a full hydrostatic pressure applied in the rock mass to a 
depth equal to the height of the dam, then, up to a certain water 
level in the reservoir, the dam will rotate in an anti-clockwise 
manner. This feature is more evident in the following case study 
of Alpe Gera dam. During the impounding of the reservoir the 
point of application of the resultant of the water load will become 
high enough so that the effect of the water pressure on the structure 
will be dominant and will then make it rotate in a clockwise direction. 

If now, we assume that we do not have any crack, i.e. that 
no full hydrostatic pressure is transmitted in depth into “he rock 
foundation, we see by looking at Figures 3.3.2 and 3.3.6 that the 
dam undergoes almost no rotation. 


In Fig. 3.3.8 are presented the results of the measurements 
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FIG. 3.3.8 DIFFERENCE IN COLUMN ELEVATIONS MEASURED WITH A HYDROSTATIC 
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(FROM BOCHKIN ET AL., 1971) 
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of the difference in column elevations for different sections of 
the dam. It is interesting to note that as soon as filling of the 
reservoir started, a clockwise rotation was measured. Having this 
in mind, we can imply that the conditions represented by Fig. 3.3.7 
constitute an upper bound in the prediction of the deformations of 
Krasnoyarsk dam. In fact, to have a clockwise rotation of the 
structure as impounding increases, full hydrostatic pressure in the 
foundation rock could only act along a small fraction of the length 
of the grout curtain. 

Using the results of the stresses computed in the gravity 
plus water load analysis and assuming a i = ] condition in the 
rock mass prior to loading, we have calculated the variation of 
the total horizontal stress, with depth, below the heel of the 
dam. The distribution of this stress together with the water 
pressure distribution, U, assumed in the calculation are shown on 
Fig. 3.3.9. This figure also shows the result of the subtraction 
of the water pressure from the total stress, i.e. the effective 
horizontal stress, oy. The region of positive stress is the 
tensile zone; the negative stresses are compressive. The depth on 
which the above stresses are acting is equal to H, the height of 
the dam. 

If we assume that the rock mass has no tensile strength, 
it would then fail down to about 50% of the length of the grout 
curtain. Therefore, we see that, in the case of Krasnoyarsk dam, 
the application of full hydrostatic pressure down to H into the 
rock mass is far too conservative. 


Consequently, it is suggested that the behaviour of 






¥o anotioee tnovst*th yo? enottsvels amfoo at eonsysttib ed? Yo 
oft To pnt lit? 26 moo2 26 jst ston of pnitasvsant ef II meh ont 
2hdy ontvsH .bevwessm 2ew notistoy setwioofs 6 ,badvese Wroviezst 





V.£.€ .pF? yd bsdneesiqey enohitbnos sit sede ylamt nso ow ,batm at . 7 | 

to enotisaroteb odd to nottotbsrq odd nt bmiod yeqqu ns sduttJenoa . 
ait to nottstoy setwioofs 6 evan ot .tost nl .meb devsyonesw 7 

edt nt sweeeva ofssteovbyd [fut .zsesetsnt patbavogmt 26 eywdauv2, oa ' 


dipnsf sit to nottoeyt [lame s profs ta6 yfno biwoa Aj07 notisbruct 
.ntstwo tuo oft to . 
utivetp oft mi baguginos esezerte ot to edtugew sid patel rm ; 
ont nt notsibnes ra 6 enimu2es bos 2faylens bso! wet6w ewlq 
to noftetysy oft betsluolso sveh ow .ontbsof of yoltq eesm A207 ~ 
sig to [aan oft woted ,Asqeb Adtw . vt ,2esrs2 fetnostronl fator ont 7 
yotsw odd dttw vorntepet 2earde add Yo notiudtss2th oat a7 
no nwode sys HoOrseluotss ond nt bemdeas . Undadetate seen 
nottosytdve ont to sfueoy sdt ewore oefs svuptt 2FaT .@.€.8 . pf 
avivastie ong .9.f .22a1te fatos ond mow? sxuz2orq WStsw ant 6 art 


ond 2f geovte avisi2og to notesy ant ‘on -2eavse [sinos tron 


fo nitgqeb eft .svi2eze,qmos sys esezsrt2 avitepan ons ienoxs oltenst — ms 
to drpten oft .H ot Tkups 27 onttos 916 eseesvte svods sd fofiw ) 
smsb ant 


efitpnarte offenses on 26h 226i ve ont sand ompe2s ow tI 4 - 
ied aneate em caembas” 7 


esol ih te foam | shone bape! ih 


— 


- 





57 


10.5 
Rock Surface 
Depth 
H: height of dam 
34.0kg/cme 22.4kg/em@ 
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PlGwaog) solO TAU STRESS, WATER PRESSURE AND EFFECTIVE STRESS 
DISTRIBUTIONS BASED ON THE RESULTS OF THE GRAVITY 
PLUS WATER LOAD ANALYSIS 
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3.4 Analysis of Alpe Gera Dam 


3.4.1 Introduction 

The Alpe Gera dam is situated on the Mallero river in the 
upper Valtellina valley in the Italian Rhaetian Alps. It was built 
to control the seasonal flows of the Scerscen and Cormor torrents 
and to produce electricity. 

The dam itself is a straight concrete gravity structure 
with a triangular profile. It is 178 meters high above the lowest 
foundation point. The upstream face has a slope of 0.03 to one 
and the downstream face, a slope of 0.70 to one. At the toe of 
the dam, a concrete block was built to increase the area of contact 
with the foundation rock. A general view of the structure is pre- 
sented in Fig. 3.4.1. 

Alpe Gera dam is founded on metamorphic serpentine schist, 


an intrusive volcanic rock. 


3.4.2 Method of Investigation and Values 
Assumed for the Study 


The rock at the Alpe Gera site has been investigated by 
two static methods: hydraulic chamber test and jacking tests. The 


following is a short description of the methods used. 


i) Hydraulic Chamber Test 


This test was performed in a gallery having a length of 
5 meters and an interior diameter of 1.7 meter. Water is pumped 
into the chamber which when filled applies an all-around pressure 
to the rock mass. Usually the deformations are measured at the mid- 


section of the chamber, by recording the variations in length of 
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four different diameters. Every test allows for the application of 
five loading cycles, gradually increasing the load up to 20 to 30 
Ka/cm*. The total duration of the test is never less than three 
hours. Figure 2.2.3 shows schematically the different features 


of the hydraulic chamber test. 


ii) Jacking Test 

The jacking test also was performed in a gallery. Very 
often, as was the case for Alpe Gera, the jacks are installed in 
the same gallery used for the hydraulic chamber test. The jacks had 
a capacity of 200 tons distributing their loads on plates having a 
diameter of 250 mm. The deformations were measured close to the 
loading plates. The configuration of the jacking test is given in 
Pigs Splits 

The available data on the tests mentioned above are pre- 
sented in Table 3.4.1. We notice that all the static tests have 
been done in the right bank at elevation 1977. The modulus of 
elasticity found varies from 45,000 Kg/em® to 256,000 Kg/cm?. The 
reason for the latter value to be so high with respect to the other 
ones is probably the fact that it has been obtained from a test 
performed in the same gallery as the hydraulic chamber tests. The 
hysteretic compression of the rock mass due to the first loadings 
has led to an increase in the stiffness value upon reloading. 

In the finite element analyses of Alpe Gera dam, two 
materials have been considered: the concrete of the dam and the 
foundation rock. 

Gentile (1964) has elaborated on the different character- 


istics of the concrete used in the construction of the dam. The 


i) 
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TABLES. 4 ai 


RESULTS OF MODULI OF DEFORMATION 


MEASURED AT ALPE GERA DAM SITE 
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Type Over- Max. Modulus Direction of Remarks 
of burden Test E 5 Measurements 
Test (m) Load 9 Kg/cm at 
Ka/cm Test Site 
18 69,000 > Only F.G. 
Heals 36 45,000 K After C.G. 
60 
100 256,000 eee Test performed 
bal le in same gallery 
100 89,000 | as the H.C.T. 
40 169,000 as 
Oete 30 
60 148,000 | 
Abbreviations 
H.C.T. Hydraulic chamber test Direction of deformation 
measurement 
Jae Jacking Test wenn nee Direction of maximum 
deformation 
F Ge Padi ingeGroucing 2© Teves ren ey. Direction of minimum 
deformation 
C.G@i Consolidation Grouting 


(From Direction des Etudes et Recherches, ENEL, Rome, Italie, 
First International Congress on Rock Mechanics, V.2, Theme 8, 


Paper 20, p. 603-615) 
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concrete has a low cement content with a modulus of 300,000 Ka/cm* 
and a unit weight of 2.6 hie Poisson's ratio for this concrete 
was assumed to be 0.18. 

As was mentioned before, the choice of the assumption of a 
modulus value for the rock will not influence the value of the overall 
calculated E. This is due to the inherent proportionality between 
modulus and displacements when using a linear elastic analysis. For 
this reason the E used in the study of Alpe Gera dam was chosen as 
100,000 Kg/cm°. Poisson's ratio was assumed to be 0.2. 

The finite element mesh used in the study of Alpe Gera dam 


is shown on Fig. 3.4.3. 


3.4.3 Results of Predicted and Measured Displacements 
i) Gravity Analysis 


No data on displacements due to the weight of the dam alone 
were available. Nevertheless, a gravity analysis was performed to 
see what would be the initial configuration of the dam prior to 
the water loading. The deformation results are shown on Fig. 3.4.4. 
We can see that, based on the hypothesis that the rock mass stiff- 
ness can be represented by a modulus of 100,000 Ka/cm*, there is 
an overall anti-clockwise rotation of the dam body as was shown for 
the case of Krasnoyarsk. The rotation of the dam is caused by a 
differential settlement of the foundation rock of approximately 


9 mm. 










Sno\pt 000,008 Fo. aufubom dim as ame aut s9x90 
stevonoo etd wot otter e'nozetot n\n 3.8 t0 diptew stim 


& to nofiqnpees st to sptodo ont biota banat dha aw 2h | 
[isvavo sit to aulsv ont soneurtAr ton Iliw Jao 93 not oul ev an Tubom 





naswisd vit i snoftyoqetg snaysdnt aft ot oub ef 2fdT .3 botsfuolso ay 
V07 .2reyfsns attests asontl 6 pate codw 23nemeostqetb bas 2utubom oe a 
26 nazora 26w msb sYsd 9QiA to ybude ont mf bsew 3 ot mnesor ety . ht 
-$.0 sd of bemuezs 25w oftsy 2'noeehot Sna\pn 000 .00F of (Na aa 

mob 6%89 9qfA to yYbuse sat nt bsey deom $hamots Sdtart oT . | i 
7 ESB. ptt no nwode ef 7 - 

= ee i 

2inemsosi gai J269M. b others to 2tfueeh €.8.6 _ on a 
incision : 

snofs msb alt to tnntaw edt ot sub ednameosigath no stieb of 7 7 . 
Oo} bemvotysq 26w 2teviens yitverp 6 229 lofaiaven adsl teve pa : 7 
ot. yotyvq msb odd to notssywertnos fstttar ont od bivow tertw 392, | : 


$.S.E pit no nwode sve edfueoy notdsmyoteb shT poi bsof vet sw sgt * ae 
-ttit2 22am A001 ont tent eteortogyrd eng no bsesd .Jedd sez n69 SW) 


at swans “no\o4 900,00f to evlubon s yd batnazerqey sd aso 22am, % 
Yo? nworle 26W 26 yybod. msb sit to noftstoy setwioolo-tins [fsvevoe ns oi : 
& yd beeuso 27 msb ont to notistoy at leveyon2sy To 9269 ont na 


| Merete Ye 1207 aotssbawot oft to hale eae 
aa 





TO oo PO 
ESPNS INNA MAIN NZ NZNNNI 
TATA AVIA | 


NN ie 


ar AWA 


: 
MAG ARID 290 silt cnentn tine w 5 F 7 


i ’ 
- 
F 
iY 7 9 
7 = - 


- 7 TW 7 | Dy ) ry 
aay 
a 








FIG. 3.4.4 DISPLACEMENTS OF THE DAM DUE TO GRAVITY ONLY 
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ii) Gravity plus Water Load Analysis (Assumption 
of a Crack below the Heel of the Dam 

The reference for the comparison of predicted and measured 
displacements in the case of Alpe Gera dam was the rotation measured 
in the foundation gallery for the whole period of impounding up to 
the normal reservoir level. The two points involved in the rotation 
measurements are point 3c and Ic. They are situated at the ex- 
tremeties of the foundation gallery, as shown on Fig. 3.4.1. 
Figure 3.4.5 shows the reservoir level and temperature variations 
together with the levelling measurements on the dam crest and in 
the foundation gallery. The eeanA instruments were located in 
block No. 10 of the dam. The period taken into account in our 
analyses was the year 1966. This was the first year when the 
water reached its normal level. The initial condition considered 
was with the water at level 2040 m. and the final one, at a level 
Ofec |Zoums 

The forces involved in these two analyses were the weight 
of the concrete, the water forces applied on the upstream face of 
the dam*, and on the reservoir bed, water forces in a crack at 
the heel of the dam penetrating approximately to a depth equal to 
the height of the structure and finally uplift forces applied on 
the base of the dam. With respect to the last forces, a maximum 
of 25% of the head was applied immediately downstream of the grout 
curtain with a linear variation decreasing to zero at mid-base. 

The displacements of the structure, with the water level 
* The upstream face of the dam was covered by steel plates. The 


modelling of the water pressure on this face is therefore very 
close to reality. 
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FOUNDATION (A) IMPOUNDMENT (B) EXTERNAL TEMP- 
ERATURE (C) LEVELING ON DAM CREST a) SECTION 
LO BDIBSECTIONS 20.) GEG MLON So 2ee (DAL EME les 
ING IN FOUNDATION d) POINT Ic e) POINT 3c 
(FROM C. TERRACINI, 1970) 
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at 2040 m. are shown on Fig. 3.4.6(a). Under the loadings mentioned 
above, the dam has tilted in an anti-clockwise direction. The dif- 
ferential settlement in the foundation rock has increased to 34 mm. 

The last analysis performed on Alpe Gera dam implied an 
increase in the reservoir level to 2125 m. The corresponding dis- 
placements are given on Fig. 3.4.6(b). 

It is obvious that between the initial and final condition 
studied here came a time when the rotation of the structure changed 
direction. This fact was mentioned in the study of Krasnoyarsk 
dam. The whole structure has therefore rotated in a clockwise 
direction bringing the differential settlements between point 3c 
and 1c back to 9 mm. 

Consequently, increasing the water level in the reservoir 
by eighty-five meters produced a net clockwise rotation cor- 
responding to a differential settlement of 25 mm. During the period 
considered Fig. 3.4.5 indicates to us a differential settlement 
of about 7.5 mm. between point 3c and lc. Therefore, our assumption 
of 100,000 Kg/cmé for the overall modulus of the rock should be 
increased to almost 335,000 Kg/cm® to account for the measured 
rotation. 

Looking back to Fig. 3.4.5, we see that levelling on the 
dam crest between the low and high level of 1966 gave a net upward 
movement of 7 mm. The different elevations calculated for those 
two conditions result in a movement in the same direction but 
having a value of 27 mm. Consequently, a ratio of 3.75 exists 
between the predicted and measured displacements. Naturally, the 


displacements of the dam crest are also influenced by the temperature 
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variations of the concrete. For the period studied the temperature 
variations are certainly not high enough to account for a sub- 
Stantial decrease of the end result quoted above. 

During the year 1968, the water level increased by the 
Same amount as in 1966, the initial and final level being also the 
same. The results obtained by the numerical analyses would then 
remain constant. It is interesting to note that the differential 
settlement between point 3c and lc in the foundation gallery has, 
in 1968, decreased to 6 mm., therefore increasing the ratio 
Enredicted! “measured to 4.15. Would this fact mean that the stiff- 
ness of the rock foundation of Alpe Gera has increased between 
1966 and 1968? If we assume that, from 1966 the foundation 
behaved elastically and we evaluate what the differential settle- 
ment between point 3c and Ic has been, for an increase of 85 m. 


in the water level we find that: 


i) in 1966 Diff. Settlement = 7.5 mm. 
id) T21nt1 967 Diff. Settlement = 6.15 mm. 
iii) in 1968 Diff. Settlement = 6.0 mm. 


This decrease in the differential settlement over that 
period suggests that, indeed there has been a stiffening of the 
foundation rock. Naturally, the improvement of the quality of the 
foundation rock at Alpe Gera does not constitute the critical 
conditions for the analysis of foundation behaviour. But, it is 
of interest to see that the competence of the rock has increased 


over the period of three years. 


3.4.4 Comparison of Actual Behaviour with the Predicted 


The overall modulus of elasticity obtained previously, for 


a 
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the conditions represented by Figures 3.4.6(a) and (b) was based 
on the comparison of vertical movements of the dam. Figure 3.4.7 
Shows the measurements of the deflection of block No. 9. It is 
assumed here that the deflection of block No. 9 is representative 
of the one which happened in block No. 10 due to their proximity. 
The curve of interest to us, is curve 'd', showing the displace- 
ments measured in the foundation gallery, at point 3c. Between 
the low and high reservoir level of 1966, a horizontal movement 

of 5 mm. has been measured. A displacement of 50 mm. is predicted 
by our analyses. Therefore, the ratio existing here would be equal 


6 Kg/cm? to account for the measured 


to 10, implying a modulus of 10 
horizontal movements. Obviously, such a high value for a serpentine 
SChist 1s rather wnrealist£1G.§ Since wthis yatiopnadibeen obtained 
assuming full hydrostatic pressure applied down to a depth equal 
to the height of the dam into the foundation rock, we have invest- 
igated the effect of varying the depth of cracking on the rotation 
and the horizontal displacements of the structure. Figures 3.4.8 (a) 
and 3.4.8(b) show the rotation and horizontal displacements res- 
pectively. The movements of the whole structure are also shown 
from Fig. 3.4.9 to Fig. 3.4.11 for different depths of cracking. 
Figure 3.4.8 indicates that a decrease in the depth of cracking 
induces a decrease in the rotation and the horizontal displace- 
ments of the dam. 

By careful examination of Table 3.4.1 we note that for 
the hydraulic chamber tests the direction in which the displacements 


had been minimum was the horizontal one, and that the vertical 


direction had undergone the maximum displacements. Furthermore, we 
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DEFLECTIONEORS BLOCK NOoe9 

(A) IMPOUNDMENT (B) EXTERNAL TEMPERATURE 
(C) DEFLECTIONS MEASURED (d) DEFLECTIONS 
MEASURED BY PENDULUM AT EL. 1981.25 m. 
(FROM C. TERRACINI, 1970) 
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see that the jacking tests have been performed in the horizontal 
and vertical direction, giving a higher horizontal modulus than the 
vertical. Therefore, both in situ tests agreed in indicating that 
the rock was anisotropic. 

Let us define the anisotropy factor as A.F. = E/E» EY, 
being the horizontal modulus and Es the vertical one. Since the 
jacking tests have been done in the directions of concern to us 
in the analysis of Alpe Gera dam, the moduli they have given were 
the ones used to obtain the anisotropy factor of the rock. The 
first test gave A.F. = 2.9, the second A.F. = 1.14 averaging at 
Avro. = 2f02 

If we now qo back to Fig. 3.4.8 and try to find the depth 
of cracking that would make the ratio of the overall horizontal 
modulus over the overall vertical modulus become 2.0, we find that 
this depth would have to be equal to 0.3H. Table 3.4.2 summarizes 
the above paragraph. 

As seen from the Table, an overall vertical modulus of 
244,000 Ka/cm? would be necessary to account for the measured 
rotation whereas an overall horizontal modulus of 492,000 Ka/cm® 
would predict the correct horizontal displacements. 

If we compare the above values witn the average of the 
horizontal moduli and vertical moduli obtained from the jacking 
tests, we find a ratio of 2.06 for the vertical and a ratio of 2.3 
for the horizontal. That is to say, the jacking tests have under- 
estimated by a factor varying from 2.0 to 2.3 the stiffness of the 
foundation rock of Alpe Gera. 


It is therefore suggested that the behaviour of Alpe Gera 
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dam corresponds to one for which a crack has developed down to a 
deoth of 0.3H into the foundation rock, the rock having a hor- 
izontal modulus of deformation of 492,000 Kg/em" and a vertical 


modulus of 244,000 Kg/cm*. 
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3.5 Analysis of Bhakra Dam 


Seo yleeeintroduction 

Bhakra dam is situated about 120 miles north-west of New 
Dethi in India, in the outer foothills of the Himalayas. The dam 
provides irrigation, power and flood control by blocking the waters 
of the Sutlej river. 

The structure is a straight gravity dam having a height 
of 225.6 meters and a base width of about 180 meters between the 
dam axis and the spillway apron section. The upstream face of 
the spillway section is vertical from the top at EL. 518 m. to 
EL. 395)m.. Below that point, it slopes at 1:20.35 for a height of 
26 m., then slopes at 1:0.88 and 1:2.0 to reach EL. 328 m. The 
downstream face of the dam has a Slope of 1:0.8. Figure 3.5.1 
shows the different characteristics of the spillway section. 

The main material encountered in the foundation of Bhakra 
dam is sandstone. This sandstone is interbedded with mudstone. 

In the vicinity of the dam there are three major claystone/silt- 


stone bands: 


i) the heel claystone band located about 30 meters upstream 
of the dam axis and dipping at an angle of about 70° to 75° degrees 
downstream. Its width varies from 30 to 45 meters. A concrete 
strut, 23 meters thick was extended from the dam to the clay band 


to transfer the loads to the upstream sandstone. 


ii) the middle claystone band having a width varying from 
6 to 9 meters. 


iii) the downstream claystone band situated under the last 
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95 meters of the spillway apron. This band does not affect the 
Stability of the dam but could influence the behaviour of the left 
power plant since this last one is founded partly on sandstone and 
partly on the claystone band. 

The position of these three major bands with respect to 
the dam is shown on Fig. 3.5.1. Several shear zones are also 


present in the foundation. 


3.5.2 Method of Investigation and Values 
Assumed for the Study 


The deformability of the foundation rock at the Bhakra dam 
Site was studied by jacking tests performed on the right and left 
bank of the river. The jack, of a hydraulic type, had a capacity 
of 200 long tons and distributed its load on an area of 920 cm?. 
The displacements were measured between the two faces of the 
gallery using extensometers capable of a precision of .0025 mm. 
The different features of the jacking test are presented in Fig. 
Secs 

Table 3.5.1 shows the results of the peering tests. The 
initial tests done gave relatively low values. The results of those 
tests are shown in the Table under ungrouted conditions. On the 
basis of the trial load analysis, which was the method of design 
for the Bhakra project, a maximum stress, in the rock, of 460 psi 
(3228 kg/cm?) had been calculated. In the Table under item 3 and 
4, for a stress range corresponding to the maximum predicted value, 
we find a modulus of deformation ranging from 498,000 psi (35,000 
Kg/cm?) tonics, COURDSHmO0eo 7/5 Ka/cm*) and averaging at 620,000 


psi (43,570 Kg/cm*) for the rock in a natural state. If we now 
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Item, Type of rock. 


3: 


Io. 


II. 
12. 


He, 


14. 


Good sandstone 
WUTAROWNOGC) og Dac aaoeas 


Good sandstone 
hNOMWER LS Soa ac aa ae 


Claystone ungrouted. 


Claystone ungrouted. 


Crushed claystone 
ungrouted 


Crushed claystone 
UNeETOUE Cee 


Shattered sandstone 
THLE WOMETROMNEG|, 5 co a oc 


Crushed sandstone 
UNLOULE CEN) sneer 


Crushed sandstone 
(After grouting but 
washing done with 
WEES? OMIA) soto on pace 


Crushed sandstone 
(After grouting but 
washing done with 
Site GYNGl AFUE) oe nia 
Good sandstone 


' (Before grouting).... 


Good sandstone 
(After grouting)..... 


Location. 


Left 
Diversion 
Tunnel 


Left 
Diversion 
Tunnel 


Right 
Diversion 
Tunnel 
Right 
Diversior: 
Tunnel 


Left 
Diversion 
Tunnel 


Left 
Diversion 
Tunnel 


Right 
Diversion 
Tunnel 


Right 
Diversion 
Tunnel 


Drift 
26h 


Drift 
20R 


Drift in 
El. 1428 
Grouting 
& Drainage 
Tunnel 
loading 
parallel 

to the 
bedding 
planes 


(*) Observations not reliable. 


VABLEMS 20 | 


Cycle. 


Ist 
2nd 


Ist 
2nd 
3rd 


Ist 
2nd 
3rd 


Ist 
2nd 
3rd 


Ist 
2nd 
3ird 


Ist 
2nd 
ard 
Ist 
2nd 
3rd 


Ist 
2nd 


Ist 
Ist 
2nd 
3rd 


Ist 
2nd 


3rd 
4th 


Ist 

2nd 
3rd 
4th 


BY THE JACKING TEST 
(FROM P.S. BHATNAGAR, S.R. SHAH, 1964) 


Stress 
range psi. 


279-841 
279-841 


279-841 
S25 
Oa 


279-1840 
279-1840 
279-1840 


279-838 
279-2524 
279-838 
279-2524 


279-526 
279-526 
279-526 


279-786 
279-786 
279-786 


279-526 
279-520 
279-526 
279-789 
279-789 
279-789 


419-838 
419-838 


525-1050 


E psi. 


0.788 xX 108 
0.72 X 108 


0.51 X 108 
0.498 X 108 
0.502 X 108 


0.564 xX 108 
0.681 X 108 
0.686 x 10% 


0.424 X 108 
0.397 108 
1.056 xX 108 
Mal 3S IKE)? 


O72 ae ele 


0.768 x 108 
0.746 X 108 
O47 10% 
0.549 X 108 
0.594 X 108 
0.074 x 108 
0.169 X 108 
ONS 2ae<ar 0% 
0.132 X 108 
0.224 X 108 
O37) SK Tee 
0.545 X 108 
0.508 x 108 
8.53 xX 108 
7.10 X 108 
7.92%, 10° 
Aes ae Os 
DO SK aeey 
4.84 X 108 
AGt rlOe 
CH ag oom fed 
6.0 x 108 

21.9 X 108 
20.1 X 108 
22.9 x 10% 


(CONTINUED) 


RESULTS OF MODULI OF DEFORMATION MEASURED 
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Item. 


15. 


17 


18. 


19. 


20. 


21. 
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Type of rock. 


Good sandstone 


(Before grouting).... 


Good sandstone 
(After grouting)..... 


Good sandstone 
(Before grouting).... 


Good sandstone 
(After grouting)..... 


Shattered sandstone 
(Before grouting).... 


Shattered sandstone 
(After grouting)..... 


Shattered sandstone 
(After grouting)..... 


Shattered sandstone 
(Before grouting).... 


Shattered sandstone 
(After grouting)..... 


Shattered sandstone 
(After grouting)..... 


Location. 


Drift in 
FE). 1428 
Grouting 
& Drainage 
Tunnel 
loading 
parallel 
to the 
bedding 
planes 
So 


Go 


Os 


Drift in 

i). 1548 

Grouting 

& Drainage 

Tunnel} 

loading 

parallel 

to bedding 

planes 
—do— 


Drift in 
El. 1548 
Grouting 

& Drainage 
Tunnel 
loading 


perpendicular 3rd 


to bedding 
planes 
Drift in 
El. 1548 
Grouting 
& Drainage 
Tunnel 
loading 
parallel to 
bedding 
planes 
—do— 


= do 


(*) Observations not reliable. 


Cycle. 


ist 


2nd 


Ist 


2nd 


Stress 
range psi. 
525-2100 


525-2100 


525-3150 


525-1050 
525-1050 
525-1050 
525-1050 
525-1050 
525-1050 
525-1050 
525-1050 


525-1050 


525-2100 
525-2100 
525-2100 


525-2100 
525-2100 
525-2100 
525-2100 
575-3150 


525-3150 | 


525-3150 


108 


I 
i=) 
(ee) 
x 


4.61 X 108 

x ro 
2.94 X 108 
3.06 xX 108 
2,59 % 10° 
22 Oe a1 0% 
25339 0od 10S 


2.84 Xe 10% 


1.085 x I0® 
MAA? Se Gy 
2165 XS10" 


2.515 X 108 
2. OF xa10° 
2 tix 10° 
r.964°° 10° 
1.98 x 108 
2507) 10% 
2.095 xX 10% 


lanl 
W 
W 
to 
KK KOS KG Os 
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o 
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look at items 13, 15 and 17, for stress ranges higher than the 


expected maximum we find a modulus of deformation ranging from 


© nsi (161,630 Kg/em-) to 5.14 x 10° psi (361,210 Kg/em-) 


6 


2 o0ex 41.0) 
and averaging at 3.95 x 10° psi (277,580 Kg/cm*). 

Of concern to the designers, were the low values found 
for the anticipated stress range. Investigations were then carried 
out to see what relative improvement could be achieved by grouting 
the rock. In doing so, it was found that to reach an acceptable 
increase in the competence of the rock, grouting had to be done 
after washing of a certain area with air and water. The results of 
the tests performed after grouting are also shown in Table 3.5.1. 
The moduli of deformation obtained after grouting ranged from 
1.98 x 10° psi (139,140 Ka/em*) to 8.53 x 10° psi (599,440 Kg/cm*) 
and averaged at 3.81 x 10° DSi) 267 4/45 Ka/cm*) . This increase in 
the modulus value allowed the designers to use a design modulus of 


deformation of 2.5 x 10° 


DS im 1754685 Kg/cm-). It is to be noted 
that the foundation rock at Bhakra was grouted up to about 15 
meters upstream of the heel, 20 meters downstream of the toe, and 
down to a depth of 15.2 meters in between. 

In the analysis of Bhakra dam two materials were taken 
into account, the concrete of the dam and the foundation rock. As 
measured by the constructors of the structure, the following 
characteristics of the concrete were entered in the numerical 
analysis: 

Unit Weight = 2.44 Tm/m 


Poisson's Ratio = Opel WF 


Modulus "of Elacticity 2 267/000 kg/cm“ 


{ 






ait nedy verptit 2opnen zeovte wot . VP bre * at ie ; 
mort pntorsy noi tenvotsb to eulubom 6 brit 4 municam besosqKs 7 
Cno\p OFS, 96) feq ®0F x 1.209 (Sno\pit 088, 18°) faq ®Or x 08.8 
(Sna\on 088,018} tq Sor x 8e.€ $s patbereve bns 
bavoi zoufsev wol silt svew  2rangiesb sii ot mysonoo +0 : 
betyys neadt syow anotispetiesval SpE egarte batsqtotans ant tot 
pnttuo“p yd bevsfios ad blo tngmavovgmt ovitelsy tedw sez oF tuo ' 
ofdetqacos ne daeoy of sett bauot e6w $F ,O2 pateb al .A907 ant 
anob sd oF bah pativoip .Wa071 sft to gonateanoes ait nr sesoront 
to etivest ofl .yedew brs y66 Adiw eons ntety93 & te pnivesw votts 
.[.8.€ ofdsT mt nwone o2fs ove paliueye vadts bomvetyeq et2zes snd 
mort bepney paiduorp yetts banbetdo noftsmrotsp to flubom sAT 
(3mo\p% QML, C82) Faq Bor x £2.38 of (no Xen OAT, RET) teq ®or x 82.1 
°or x (8.€ #6 bepsrevs -bns 
to eulobom apfzsb s sev of avenphesh oft bewol fs aulsv 2ufubon sft 


G ¢ 
af s269yont 2hnT .("mo\pX ght. §eS) fae 


beaten od ot af ST mo\or 280,a0f) 2g %or x &.8 to aoftenvotsb | 
ef tuocs of qu batuoyp zew syierd te. door noitabavet ont Seat 
bis .90¢ oft to mesytenwob zysd9n O8 .f 95n att 40 ‘neat equ oie | 
ngawied nb avasom S.2T to itqsb 6 od nwob - 
nsasy svaw efeinsiem owt meh sybsld to: 2taylans alt ni. : 
2A .A204 noftebnuo? add bas mb sid +o ateron09 orld es otnt | 
enitwolfot ont .oxusourte eit to evotoviteno9 al vd bowuesom hi! 


CA4he 


Isotyomun oat tt barodna, 1a seein 





88 


With respect to the rock, the same assumptions as for the 
Studies of Krasnoyarsk and Alpe Gera dam were made i.e. 
Unit Weight = 0.0 
Poisson's Ratio = 0.2 
Modulus of Elasticity = 100,000 Kg/cm@ 
The finite element mesh used in the study of Bhakra dam is 
shown in Fig. 3.5.3. The maximum spillway section, at the center of 
the dam, was the one which was analyzed. This section is section 


#20. 


3.5.3 Results of Predicted and Measured Displacements 

i) Gravity Analysis 

A gravity analysis of the structure was eideeceren to see 
what would be its configuration prior to water loading. The displace- 
ments of the dam are shown in Fig. 3.5.4. 


ii) Gravity plus Water Load Analysis (Assumption 
of a Crack below the Heel of the Dam 


The observation of the behaviour of Bhakra dam consisted of 
measuring the deflections of the structure, recorded by plumb lines, 
and to measure settlements of benchmarks fixed in a gallery at 
EL. 373.4 m. The horizontal deflections were measured relative to 
a point in a gallery at EL. 336.8 m. which was assumed to be fixed. 
Due to this fact, horizontal movements of the structure were not taken 
into account in the present study. Permanent settlement readings were 
started in December 1963 and continued ever since. Impounding had 
started in 1958 which means that up to 1963 the foundation rock has 
been subjected to cyclic loading and unloading by the reservoir water. 
Figure 3.5.5 shows the settlements of six benchmarks together with the 


variation of the reservoir water level from December 1963. In the 
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FIG. 3.5.4 DISPLACEMENTS DUE TO GRAVITY ONLY 
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FIG. 3.5.5 SETTLEMENTS OF BENCHMARKS 
(FROM BHATNAGAR ET AL., 
1967) 
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following analyses we will investigate the variation in settlement 
from low water level, EL. 443.5 m. (1455 feet) to the high water 
level of 1964 at EL. 506 m. (1660 feet). 

Up to August 1963, the dam behaved normally. In the period 
including September and October 1963, extensive leakage was measured 
through cracks in the transverse strut gallery. Two pipes had 
been installed for uplift pressure measurement at the concrete-rock 
contact at the upstream extremity of the strut. These two pipes 
started registering full hydrostatic pressure on August 30, 1963. 
The observation of cracks in the transverse strut gallery and the 
leakage indicated that separation of the strut from the dam had 
occurred at or near the longitudinal joint plane, 23 meters up- 
stream of the dam axis, in a direction normal to the strut surface 
Slope. In our analyses, since the strut was separated completely 
from the dam body, a crack was assumed through it where full hydro- 
static pressure was developed. Also, full reservoir head was applied 
around the strut, on the rock, up to approximately the heel of the 
dam. Figure 3.5.6 shows the loading conditions on the structure. 

As mentioned above, the dam was loaded with the low and high 
reservoir level of 1964. The water reached elevation 506 m. (1660 
feet) at the end of September 1964. Table 3.5.2 shows the measure- 
ments of uplift pressures, expressed as a percentage of net head, 
for different sections of the dam. As we can see, the uplift at 
section #20 was practically non-existent. For this reason, uplift 
loading on the base of the dam was neglected in the analyses. 

At this stage, consideration has to be given to the fact 


that the analyses assumed plane strain conditions even if Bhakra dam 
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FIG. 3.5.6 LOADING CONDITIONS ON THE DAM AND FOUNDATION AFTER 
SEPARATION OF THE STRUT 
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is located in a narrow gorge. At section #20, a deflection parallel 
to the dam axis of the order of 2.5 mm. was measured. This displace- 
ment was directed towards the left bank. It is not felt that such a 
small displacement is enough to abandon the plane strain assumption. 

By looking at the mode of vertical displacement shown on 
Fig. 3.5.5, it was expected that, at the most, a shallow crack, below 
the heel of the dam could develop. We see that for the water level 
variation to which the structure was subjected, settlements increase 
as the water level increases and vice-versa. In the past experience, 
i.e. the studies of Krasnoyarsk and Alpe Gera dams, we saw that, for 
a water level range going from about 50% to 90% of the dam height, 
settlements changed direction. This was based on the assumption of 
a crack in the foundation rock, below the heel of the dam. Bhakra 
dam departs from this behaviour. Cracking was therefore assumed 
down to a depth of 0.35H as a starting assumption. Figure 3.5.7(a) 
and (b) show the settlement of the structure under the specified 
loading. At low level, the displacement corresponding to bench- 
marks B-7 and B-8, situated in the plane of the dam axis, was 135 mm. 
Raising the reservoir level gave us a displacement of 131 mm. 
Therefore the dam has moved up by 4 mm. 

The second analysis assumed a crack down to 0.15H in the 
foundation. The displacements are shown on Fig. 3.5.8(a) and (b). 
The low level gave a vertical displacement of 131 mm. and the high 
level a displacement of 136 mm. On the basis of this assumption 
we had an increase of 5 mm. in the settlement value. 

The final analysis looked at the conditions where no crack 


was present in the foundation rock. The displacements corresponding 
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to this situation are shown on Fig. 3.5.9(a) and (b). Under the 
low reservoir level, a vertical displacement of 128 mm. was pre- 
dicted. Increasing the reservoir level to EL 506 m. gave a final 
settlement of 142 mm. Under the conditions of no crack, we find 
that the dam is subjected to an increase in settlement of 14 mm. 
For benchmarks B-7 and B-8 which were situated in the 
middle of the dam, Fig. 3.5.5 gives us an increase of 4 mm. in 
the vertical displacement. Comparing this value to the variation 
of settlement obtained for the second analysis (Depth of Crack = 
0.15H) we find a ratio of 1.25 between the predicted and the 
measured. If we now compare the measured settlement to the one 
obtained in the final analysis (No Crack), we find a ratio of 3.5 
between the predicted and the measured. The ratios of 1.25 and 
3.5 give an overall value for the modulus of elasticity of 125,000 
Kg/cm® and 350,000 Ka/cm® respectively. Comparison has now to be 
made with the in situ measured values. In Section 3.5.2 we have 
quoted a value of 43,570 Ka/ em? for the modulus of deformation of 
the rock in a natural state. Consequently E 


overall/ "measured : 


2.87 for a depth of cracking equal to 0.15H and E s/E 


overall’ measured 


8.0 for the no-crack analysis. 
Comparison of the same moduli with the ones measured after 


grouting and averaging at 267,745 Kg/em® gives the following ratios: 


E = 0.47 for crack = 0.15H 
overal] 
measured 
and 
EveraT| = 1.30 for the no-crack analysis 
Eneasured 


Table 3.5.3 summarizes the above results. 
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caw Comparison of Actual Behaviour with the Predicted 


As we have underlined in Section 3.5.2 extensive con- 
solidation grouting was done in the foundation of Bhakra dam. 
Measurements of uplift pressure on October 3, 1964 showed that no 
uplift had been recorded in the foundation rock, downstream of 
the heel. It is believed that the quality on the consolidation 
blanket was high enough to prevent cracking and subsequent prop- 
agation of uplift pressure in the rock. 

Based on these facts, it is suggested that the conditions 
represented by the final analysis, in which the presence of a crack 
was ignored, are more likely to represent the actual behaviour of 


Bhakra dam. The E based on an average E 


overall "measured Fall, 
after grouting, therefore gives us a value of 1.30. 


It is of interest to discuss the value of E which 


grouted 
is used here for comparison purposes. Grouting was carried down 


15.2 meters into the foundation rock. The change in stress caused 

by the dam itself and the water loading is responsible for the 
deformation of the rock. This change in stress does not affect 

only a layer of 15 meters but is felt down to great depths. The 
modulus of deformation of the rock immediately below the consolidation 


blanket is likely to be smaller than E Mat 1s to say, cnac 


grouted. 
the in situ E we compare our overall modulus obtained by the 


numerical analysis to should include the composite effect of a stiff 
layer and of a region of expected lower modulus. This composite E 


would be smaller than E but higher than E There- 


grouted ungrouted. 


ratio 


fore, it would be legitimate to ah an EO VErATT cin cath 


higher than 1.30. 
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For the Bhakra project a modulus of deformation of 175,685 
Kg/cm® was used for desiqn. We then see that the decision of the 
designers brought a factor of safety of 2.0 on the deformability of 
the rock mass at Bhakra dam site. 
3.5.5 Influence of Poisson's Ratio 
on the Final Result 

In the available data on Bhakra dam, values of Poisson's ratio 
for the different materials encountered in the foundation had been 
measured. The results of those measurements are shown in Table 3.5.4. 
In the petrographic classification within the table, calcareous 
subgreywacke corresponds to what was called sandstone in the field. 
Since this material was, as mentioned before, the main material under 
the dam, it is the one we will refer to in the following. 

As we can see in Table 3.5.4, Poisson's ratio was measured 
for air dry and moist conditions, for different stress ranges going 
from 0 to 200 psi (0 - 14 kg/cm’) to 0 to 1000 psi (0 - 70 Kg/cm*) 
for the air dry conditions and from 0 to 500 psi (0 - 35 Kg/cm*) to 
0 to 4000 psi (0 - 280 Kg/cm®) for the moist conditions. For the 
tests performed under air dry conditions, Poisson's ratio takes a 
minimum and a maximum value of 0.0 and 0.17 respectively. Under 
moist conditions the minimum Poisson's ratio found was 0.003 and 
the maximum, 0.13. 

Obviously, rock under a dam will have a water content 
greater than zero. The values given for a 75% moisture condition 
are therefore more realistic with respect to the field conditions. 

In the finite element analysis, the effect of varying 


Poisson's ratio, on the final result of the previous section, was 
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studied. A Poisson's ratio of 0.1 and 0.4 were input to the analysis. 
Table 3.5.5 summarizes the results obtained. Figure 3.5.10(a) shows 
the variation of the overall modulus of deformation needed to account 
for the measured displacements as a function of Poisson's ratio. 


Fig. 3.5.10(b) shows the influence of Poisson's ratio on the Ee yera he 


E ratio. 


measured 


For the range of values measured on samples taken from the 
field, we see that there is no appreciable difference in the overal] 


modulus and consequently in the E Yate Tacit, 


overall’ ‘measured 
with the assumption of a Poisson's ratio equal to 0.2 in our analyses, 


we are on the conservative side. 
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FIG. 3.5.10 INFLUENCE OF POISSON'S RATIO (a) ON THE 
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3.6 Analysis of Underground Powerhouses 


The only powerhouse that will be analysed is the one excavated 
at the Oroville dam site, in Northern California. This underground 
opening has already been studied by Professors J. M. Raphael and 
R. W. Clough at the University of California. Results of its behaviour 
have been presented by Kruse (1971). One reason for redoing the 
analysis of the Oroville powerplant is to see if the procedure which 


will be used here will reproduce the results already obtained. 


3.6.1 Boundary Conditions 


The boundary conditions used ‘in the analysis of Oroville power- 
plant are as shown on Fig. 3.6.1(a) and (b). Further explanation of 
those boundary conditions are given in a subsequent section. 

The top of the mesh was sloping, modeling the flank of a 
mountain. Its average height and width were about ten times the 
width of the underground opening. 

On the vertical sides of the mesh movements were allowed in 
any direction. The bottom boundary of the mesh was completely fixed. 

Simulation of the excavation, which constitutes the loading of 


the finite element mesh, will be described in the following section. 


3.6.2 Procedure Used in an /inalysis 

In the analysis of are tndererotnd opening, knowledge of the 
initial state of stress in the rock mass is of primary importance. 
If from its determination it is found that the horizontal stress 
is related to the vertical by Poisson's ratio, a gravity analysis, 
in which the horizontal displacements of the side boundaries are 


fixed, will give the measured state of stress. On the other hand, 
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if the above relation does not exist, it is probably due to local 
heterogeneitiesin the rock mass or to the tectonic forces which 
depend on the earth's crust activity. 

Therefore, knowledge of the initial state of stress allows 
us to load the lateral vertical sides of the mesh in order to generate 
the stresses in the finite element model. Naturally, if this is 
done, movements of the lateral boundaries must be allowed in any 
direction. 

In the study of the Oroville powerplant two analyses will be 
performed. As explained above, the first one consists of inducing 
the initial state of stress in the model. In the same analysis, 
the forces which will simulate the excavation of the rock mass can 
be calculated by applying the following finite element equation, 
(Zienkiewicz, 1971): 

{rhe - J Cel des) dy (1) 
Where: iF = is the matrix corresponding to the nodal 
forces for one element, 
[eB] = is the matrix representing the geometry 
of an element, 
(ooh = is the matrix corresponding to the initial 
stresses in one element, 


is the volume of one element 


= 
il 


If the calculation involved in this equation is done for all 
the elements inside the outline of the opening and having a side in 
common with this outline, nodal forces will be calculated all around 


the opening. 
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It is important to note that the state of deformation the 
finite element model is in, at the end of this analysis, is assumed 
to represent the state of equilibrium in the rock mass, before starting 
the excavation. 
In the second analysis the following items are applied: 
i) the same boundary conditions all-around the mesh are kept 
but no side forces are applied, 
ii) for all the elements included inside the outline of the 
opening, properties of air (E*o) are given, 
iii) the weight of the material surrounding the opening is set 
equal to zero, 
iv) the forces found in ine first analysis are applied to their 
respective nodes, around the outline of the opening. 
Performance of this second analysis will give the displacements 
caused by the excavation itself. The analysis being a linear elastic 
one, the final state of stress around the opening is obtained by 


adding the states of stress of the two analyses. 
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3.7. Analysis of Oroville Underground Powerhouse 
3.7.1 Introduction 


Oroville underground powerplant is situated in the state of 
California, in U. S. A., about 85 miles north of the City of Sacramento. 
Oroville dam and powerplant are part of the State Water Project in 
California, which has the objective of distributing water from 
northern regions to the southern regions. 

The powerplant itself contains six generating units. Its 
length is approximately 170 m. (550 feet). Its maximum width and 
height are 21 meters (70 feet) and 42 meters (140 feet) respectively. 

The underground chamber has been excavated in metamorphic rock 
high in amphibole minerals. The rock was called amphibolite. At 
the excavation site three joint sets charaterize the rock mass. | 
Shears and schistoze zones are Neots ene A more detailed descrip- 
tion of the site geology is given in Kruse (1971). 

Excavation of the opening started in March 1964 and was 90% 


completed by May 1965. 


3.7.2 Methods of Investigation Used and Values 
Assumed for the Study 


i) Deformability Characteristics 
A publication by Kruse (1969) discusses the different methods 


used in the determination of the modulus of deformation of the rock 

at Oroville dam site. Table 3.7.1 summarizes those results. As can 
be seen a very extensive testing program was undertaken. The results 
we are concerned with here are those of the flatjack tests, the 

plate bearing tests (or plate jacking tests) and the tunnel relaxation 


method. 
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The primary objective of the flatjack tests performed was the 
determination of the initial state of stress in the rock mass. How- 
ever, the flatjack test can also give an estimate of the modulus of 
deformation knowing the stress transmitted and the displacement under- 
gone by the grouted pins. Thirty measurements with the flatjack 
test gave moduli ranging from 100,000 Kg/em® (164ex 10° psi) to 
1,150,000 Kg/cm* (16.4 x 10° psi) and averaging at 530,000 Kg/cmé 
baSex 10° psi). 

Five plate jacking tests were done. Their orientation being 
mainly perpendicular to the surface of the rock in a gallery, they 
are likely to be more influenced by the destressed zone around the 
Opening. Low values of the modulus of deformation would then be 
expected. Table 3.7.1 gives values ranging from 85,000 Ka/cmé 


(1.2 x 10° psi) to 125,000 Kg/cm* (1.8 x 10° 


psi) and averaging at 
105,000 Kg/em* (1.5 x 10° psi). 

The tunnel relaxation method consists of computing analytically 
the displacements around an opening and to compare these displacements 
to the actual measured ones. The position of the measuring device 
being known, with respect to the opening, displacements can be cal- 
culated and plotted for a range of values of the modulus of deform- 
ation. Upon excavating, the measured displacements are also plotted 
on the same graph and a best-fit line can give the modulus of deform- 
ation based on the measured disp\acements. This method has the 
advantage of considering a very large volume of rock in the comparison. 
However, it is limited with respect to the shape of opening that it 
has been derived for. More details are given in Kruse (1969). 


Twenty-two tunnel relaxation comparisons were made and gave 
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moduli ranging from 42,000 Kg/cem* (0.6 x 10° psi) to 530,000 Kg/cmé 
(7.5 x 10° psi). The average modulus of deformation was 180,000 
Kg/em* (2.6 x 10° psi). 

From the above data, it is seen that a considerable scatter 
exists in the modulus value. However, since the tunnel relaxation 
method and the plate jacking test load large masses of rock, their 
results were thought to be more representative. For the Oroville 
underground chamber a design modulus of 105,000 Ka/cm® Glauy 10° 
psi) had been used in the previous analyses. The same value was used 
in our analysis. 

No data was available for the value of Poisson's ratio of 
the rock at the Oroville site. It was assumed equal to 0.25. 

The weight of the material was taken equal to 2.9 Tm/m>, 

11) InsSituelnitialsStresses 

The determination of in situ initial stresses at the site 
of the powerhouse was done by two methods (Merril et al., 1964): 
the flatjack and borehole deformation methods. 

The flatjack method performed at the Oroville site was 
similar to the one already described. The reference for displace- 
ment measurements was an assembly of eight pins, four grouted on 
each side of the flatjack. 

The borehole deformation test used an overcoring technique. 
Measurement of deformations upon overcoring was made by a gage 
inserted in a drilled hole. The deformation consisted of a change in 
diameter at a specific orientation. Therefore consecutive measure- 
ments at three different depths from the rock surface needed to be 


done in order to obtain a complete state of stress. 
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In spite of the disadvantages of each method, their agree- 
ment was relatively good. From the measurements it was concluded 
that the stress field was essentially hydrostatic and equal to 
35 Ka/cm* (500 psi). This stress corresponded quite well with the 
weight of the overburden. In our analysis the horizontal stress 
was then made equal to the vertical stress determined from the 


weight of the overburden. 


3.7.3 Results of Predicted and Measured Displacements 


In the analyses performed on the Oroville powerplant, only 
the rock overlying the opening has been considered. The powerplant 
was built under the upstream shoulder of the dam. As was mentioned 
earlier, the excavation was almost completed in May 1965. At that 
time there was not a considerable thickness of material directly 
on top of the chamber. The finite element mesh used in the study of 
Oroville powerplant is depicted in Fig. 3.7.1. 

Complete results of the deformations measured around the 
chamber are given by Kruse (1969) and Kruse (1971). Extensometers 
had been installed in the arch and the walls of the opening. They 
were approximately 6m (20 feet) long in the arch and 12m. (40 feet) 
long in the walls. Figure 3.7.2 shows their relative position. 

In the arch of the opening, thirty-nine deformation measurements 
were made ranging from .05 mm. (.002 inches) to 1.5 mm. (.059 inches) 
and averaging at 0.30 mm. (.012 inches). For a depth of six meters 
from the rock surface our analysis gives a relative displacement of 
1.33 mm. (0.52 inches). Therefore to account for the measured dis- 
placement of the arch the modulus used should be multiplied by a 


factor of 4.45. This would result in a modulus of deformation equal 





ees cane 
ar a oi 


~aetps viadd <bodktem toss: 46 saldehsiaeeaaialiagaieal ' ™ 
bobufonos 2aw tf etmemanyessn, od? mova, somatennieintaentt 
od [pups bws attetzovbyd yi iahiaeees esw dfott 2eatte ond tai 

and dttw Ifew stiup bebmoqzeyves 2eowta etdT .(feq O0e) Sino \ pd a 
2eante fetnostvod ond afeylisns suo A imabyudisve sit to tdptow 
903 mov? bentrystsb ezesi2 Tsabivev ent od feups abam neds 260 
nabiediave ent to tdptew ‘ 


atoamooa teat bowaseM bos bedotband Yo aftuesd €.%.€ 








vino ,inslavewor offfvorv0 srid wo bemvoTisg zs2y[6n6 odd al 
jnsfarswoq sat .bsysbteneo nsed 2aii ont nego ant ontyfrave lox sat 
bonotinem 26w 2A .meb sett to ybhivode meaitequ ott yaonu sftud 26w 
jeg SA .2d8F ysM nf badsiqmoo Feomis eew not t6vs2x9 ont at i163 
yisasiib fsivetem Yo 2esndotng sldsvebrenos & ton 26W syed omits 
to vbut2 ot nt beeu dzem tnomefo ettnt? sAT .qodmerds Sat to got no A 
.[.N.€ .pFF af betofaeb ef Jnsfarewoq i 
ond baveis bewwesem enotssarvoteb. sit to 2ifuesy stoiqmod 
2retemoenatx3 .( Tl) sauwt bos (R88T) seuvi ‘vd NOVED S16 ~vedmarts 
yaiT .painaqo oft Yo ellew add brs nove sd at bof stent nesd ber 
(tot O) .mS{ brs doves ont at onol (Soot 0S} ma ylstemtxorqas S196 
Wortteoq svitstsy yredt aware S.3.8 otupt4 -alisw aft nt prot 
2dnamewwesam Not smoteb ontaqyirtis -ovtnago silt to rove oft nl ran 


eure a as ne geet cpaplant. 













NN 
hwy 
Sy 
SN 
Z 





MA 


SAAAVYVA 
Z 


AZ 
KD 


INNES 
/\ 


SSREe 
WV. 


RAW 


ns 


» 


a 


Oe 


WiVAVAV A: 






NS 
VI \VVININNNXSS 
NIN 
oS 
N 


S 
— 
Za 


We 


VV 
a 
te 

Ara & 


KS 

oe ea es SESE 
ages 
ZS ESE 


\ 





ZEN 


WA 
N 
VV\ 


= \A\ 
VAZAVAVA 


WN 


Lv 
Se 


Wi 


SS 


y 
WIAKKAARNSSS 










VY U TTIKSS 

Ly Ye ZAKS 
VLAN 
UY AKNN 







MESH BUSE 
UNDERGROUND POWERPLANT. 


ly 





FIG. 3.7.2 RELATIVE POSITION OF EXTENSOMETERS FOR WHICH 
DISPLACEMENTS ARE QUOTED. 
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to 470,000 Kg/cm* (6.5 x 10° psi) for the arch portion. The finite 
element analysis performed at the University of California gave a 
relative displacement of 1.65 mm. for the arch of the opening. This 


last displacement would result in a modulus equal to 580,000 Ka/cmé 


(Gr2ex 10° psi) for the arch portion. 

At the end of the excavation the displacements measured by 
the horizontal extensometers ranged from 1.38 mm. (0.054 inches) to 
6.3 mm. (0.247 inches). The average was 2.9 mm. (0.114 inches.) 
The relative displacement for a distance of 12 meters and calculated 
by the numerical analysis was the same for the right and left wall. 
Tt was equal to 2.7 mm. Almost no difference exists between the 
calculated and measured displacement. For all practical purposes 
it can be said that a modulus of 105,000 Kg/cm* (1.5 x 10° psi) 
correctly represents the walls'deformation behaviour. The relative 
displacement calculated at the University of California was equal 
to 3.55 mm. (0.140 inches). Although resulting in a higher modulus 
than the one we have calculated, this relative displacement is further 
from the actual measured one. All of the above results are summarized 
in Table 3.7.2. The absolute deformation of the rock surface, cal- 
culated in our analysis, is also given in Fig. 3.7.3 

The fact that such a big difference exists between the moduli 
of deformation needed to account for the measured arch and wall move- 
ments is of particular importance. Some possible reasons, were sug- 
gested in Kruse (1971). The one which was considered as the most 
plausible was the following: "The difference in bolt spacing and 


installation specification may account for the apparent difference 


in crown modulus and wall modulus." 
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FIG. 3.7.3 ABSOLUTE DISPLACEMENTS OF THE OPENING. 
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Indeed the rock bolts of the arch were installed in a four feet 
center to center pattern whereas the bolts of the walls were six feet 
away from each other. 

To summarize, our analysis of Oroville powerplant grossly 
agrees with the one performed at the University of California. 


However, our results show a closer agreement to the actual behaviour. 
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CHAPTER IV 
CONCLUDING REMARKS 


Orne Deformability of Rock Masses 


In this study we have tried to prove that the overall stiff- 
ness of a rock mass subjected to the loading of a dam and adjacent 
reservoir, for example, is higher than the stiffness arrived at by 
an in situ static test. It is important to note that for the three 
cases analysed, different moduli (modulus of deformation or modulus 
of elasticity) have been used to obtain the final ratios. In the 
case of Krasnoyarsk dam the in situ modulus was found by subjecting 
a plate to vertical loading. The modulus obtained from such an in 
Situ test could be regarded as a deformation modulus. In the 
analysis of Krasnoyarsk dam the measured displacements used for 
comparison purposes were given with respect to the beginning of the 
construction of the dam. Therefore these displacements included the 
effect of gravity, water loading, creep of the rock mass etc. The 
analytical overall modulus obtained after comparison is therefore 
an overall modulus of deformation. Care has to be taken so that the 
overall analytical modulus obtained and the one determined from the 
in situ test are of the same type. For Krasnoyarsk dam an overall 
modulus of deformation was then compared to an in situ modulus of 
deformation. 

The study of Krasnoyarsk dam led us to the conclusion that 


an overall modulus of deformation of 210,000 Ka/eme is representative 
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of the measured structure behaviour. The average in situ modulus 
of deformation was equal to 100,000 Ka/em@. 

The case of Alpe Gera dam was different than Krasnoyarsk 
dam. Here, the displacements used for the comparison had been measured 
long after the first impounding of the reservoir. The structure had 
been subjected to cyclic loading. Therefore the analytical overall 
modulus obtained was an overall modulus of elasticity. Consequently 
this modulus has been compared to an in situ modulus of elasticity. 

In the analysis of Alpe Gera dam, we found that the overall 
moduli of elasticity of 492,000 Ka/cm* and 244,000 Ka/em? would 
account for the measured horizontal and vertical displacements res- 
pectively. These moduli average at 368,000 Ka/em”. The hydraulic 
chamber test performed at the dam site gave an in situ modulus of 
69,000 Ka/cm*. The plate jacking test taken to a stress level of 
100 Ka/em* averaged at 175,000 Keene 

Finally, the case of Bhakra dam is similar to Alpe Gera dam 
in the sense that the structure has been subjected to many cycles 
of impounding. The analytical overall modulus, with respect to our 
terminology would then be an overall modulus of elasticity. In the 
section on Bhakra dam the moduli determined by the static tests were 
called moduli of deformation. The overall modulus of elasticity, 
obtained from the analysis, was compared to a modulus measured after 
grouting had been performed. The very high values obtained after 
grouting suggest that the modulus of elasticity of a test would be 
very close, if not equal, to the quoted modulus of deformation. 

In the study of Bhakra dam an overall modulus of elasticity 


of 350,000 Ka/cm* was found to account for the measured displacements. 
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The average of the in situ tests moduli, measured after grouting, 
was 267,745 Kg/cm. 

In this work only one underground opening was analysed. The 
Study of Oroville underground powerhouse revealed that moduli of 
470,000 Kg/cm® and 105,000 Kg/cm? were necessary to account for dis- 
placements of the arch and walls respectively. The average of the 
moduli determined by the plate jacking test is 105,000 Kg/cm°. 

Figure 4.1.1 summarizes the results obtained all through this 
study. Ine data of Fig. 2.3.2 1S also* included. 

From Fig. 4.1.1 it can be concluded that large in situ static 
tests will underestimate the value of the modulus of deformation or 
elasticity of rock masses. This can readily be understood: 

If we refer firstly to the case of a dam, we find that large 
in situ static tests are usually performed at shallow depths from the 
surface. At that depth weathering has a more marked influence on the 
rock; very often surface excavation of rock is undertaken which 
results in a stress release and subsequent opening or widening of 
joints; in areas of severe climatic conditions, freezing of water 
in the winter plays a role in increasing the deformability of the 
upper strata; in situ static tests being usually performed in gal- 
leries, the surrounding destressed zone will affect considerably the 
value of the modulus obtained from the test. On the other hand, 

a structure like a dam will carry its load down to great depths in 
the rock mass, where the rock is usually sounder. 

In the case of an underground opening all of the above factors 
might affect the value of the modulus determined from a test. But, 


since underground openings are often excavated deep within the rock 
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body, the last factor mentioned above is the one that plays the 
dominant role. Therefore, if provision is made for long buried gaaes, 
to go beyond the zone of disturbance in a test adit, it is likely 

that the modulus obtained from a static test will be more represent- 
ative of the rock mass. 

From the studies of many dam sites Snow (1968) has evaluated, 
from permeability measurements in drill holes, the fracture spacing, 
the size of fracture opening and fracture porosity of rock masses. 

His conclusions are: 

i) the fracture porosity decreases logarithmically 
with depth, 

ii) fracture spacing increase with depth, 

iii) fracture openings decrease with depth, 

iv) neither spacings nor openings are notably different 
from one rock type to another, nor are porosities or permeabilities 
which depend on the first two. 

All of these items are consistent with an increase of rigidity 
of the rock with depth. 

Rocha (1969) has measured with the LNEC dilatometer an increase 
of modulus with depth. His results are shown in Table 4.1.1. 

Burland and Lord (1969) have measured, under ideal conditions, 
an increase of the modulus of elasticity with depth. Their results are 
SHOWN Neds 4012.2. 

Many other investigators (Hast, 1967; Voight, 1969; Li, 1970) 
have measured an increase in the horizontal stress with depth. Figure 
4.1.3 shows some results gathered by Voight (1969). This increase in 


horizontal confinement with depth would also result in an increase of 


a 0) oe an 






ee et bod 
. 29060 bstvud pitof 4Ot sham 2h noteiverq: +t sroteratT ia ‘anima 
volt! 2i dF bs Seat 6 nt sonsdyusetb to anos silt bao ad 09 of 
-Insesyqoy atom od liw teed sttete 6 mov? banteddo 2uTubon ons soc y 
| -2esm Aooy ant to evtts 
badsulsve zed (88CT) wone zett2 mab ynem to esfbuse odd mot? | ms 
-bnifosge sytosnt oft -2afod [fivh at 2dnemewesem yt tl tdsomieq mov? P 
.29226m ASON To yITzoy%og syuTae Xt bas enhveqo swusostt To esta orld 
:ovs enoteufonoo ett 


6® 


Vi leotard trspol 2azso10ab yst2ovog Studostt ont (t 


ettqab Ad tw 
-fqeb Atiw sessront ontosue opines? (ti 
ef3qeab dtiw sessyosb epnfnaqe ous 2617 (tit n 
tnowtiib yfdsson 916 2ontnaqo ion zhiroege rend ton (vt . . 


esitiltdsemraq 10 esrifeovg o16 von ,vedtons of sqyt i201 ano aot 
owt teat? sad ao bneqab fofrw 
\“itbipty to sessvont ns dttw tastetencs ss amet e2ens Yo TIA : 
-Haqeb ddiw dso1 st to 
s2ssrori ns vetemotsifb ISNJ odd Adtw bowwesam zor (88@T) sitoos Par 
-F.1.8 sideT mt nwote sis 2ttueoy ety -dgab: cist 2ufubon to 
eenorsrbnos fesbi wbay ,bowesem aven (edef) bro) bas base 
916 2tiues1 aiedT .Atqeb natw vttottests to euTubom end. to 9289720 ns 
Shee “Br nt mwore 
bial td ‘aon tapton BOE «3eni8) nena ce. ¢ 


ry, _Atogh Ht > ; 
oe ae Gee Ot ae io : 
7 : View - ahs 
wen > 2 Sat p x i 
p Sar ort : T\ 
E * om 


7 
ee aD 






TABLE 4.1.1 


VARIATION OF THE MODULUS OF 
WITH DEPTH 


(FROM ROCHA, 1969) 


DEFORMATION 


Modulus of 
Deformability 
(kg/cm@) 
176,500 
1555550 
181,250 
201,500 
2325/50 
258,250 
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COMPARISON OF YOUNG'S MODULUS E DERIVED 
FROM THE PLATE TESTS AND LABORATORY 
TESTS WITH VALUES DERIVED FROM FULL- 
SCALE TANK TEST 


(FROM BURLAND AND LORD, 1969) 
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FIG. 4.1.3 AVERAGE HORIZONTAL STRESS AS A FUNCTION OF DEPTH 


DATA FROM SEVERAL SOURCES (CF., HAST, 1967; 
VOIGHT, 1967, 1969). VALIDITY OF SOME POINTS 
MAY BE QUESTIONED (E.G. JAEGER, AND COOK, 1969 
FlGee 14% 3) 


(FROM VOIGHT, 1969) 
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modulus. 

All of the above mentioned references clarify the fact that a 
higher modulus than the one measured by an in situ static test is 
needed to account for the measured deformations of a structure. 

It will certainly take many more analyses like the ones which 
have been performed here, to be able to set a criteria regarding the 
overall deformability of rock masses. But, at this stage, more 
confidence may exist in the mind of the designers of large structures 
knowing that the overall stiffness is likely to be better than that 


based upon common in situ tests. 


4.2 The Behaviour of Structures 

Since this work has dealt mostly with concrete dams, the 
following observations are of interest. 

It has been recognised by many (Zienkiewicz, Toeseevernans 
1966; Serafim et al., 1967; Sabarly, 1968; Tizdel, 1970) that the 
state of stress induced in rock foundations upon impounding is of 
primary importance in the evaluation of the subsequent deformation 
behaviour of the structure. As is mentioned in the references above, 
a tensile state of stress will exist at the heel of the dam when its 
reservoir is filled. This tensile state of stress will make the 
natural fractures of the rock open and let full hydrostatic pressure 
enter the foundation. Depending mostly on the pre-existing stresses 
in the rock mass, a crack will propagate until these initial stresses 
counteract the effect of the water pressure. 

Very recently the barrage of Malpasset in France has failed. 
In a report presented by the Commission charged to study the catastrophe 


(1966), it was concluded that tensile failure of the rock upstream of 
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the dam and subsequent propagation of full hydrostatic pressure in 
the rock mass were responsible for the initiation of the failure 
of the dam. 

Although of extreme importance, with respect to the behaviour 
of a structure, seepage forces and full hydrostatic pressure in a 
crack at the heel of a dam have been rarely considered in its design. 
Rocha (1970) has analysed the effect of seepage pressure in the 
foundation rock below a concrete dam. He has used the finite element 
method and has also done model studies. His results show that water 
pressure applied on the grout curtain of a dam has a significant 
effect on the stresses and deformations of the structure. 

Throughout this study, the concept of introducing full hydro- 
Static pressure down to a certain depth, at the heel of a dam, has 
been explored. For the three dams analysed, it has helped in under- 
standing their behaviour. 

At Malpasset the conditions led to a catastrophe. One of the 
recommendations of the Commission was to assure that propagation of 
full hydrostatic pressure at the heel of a dam would not take place. 
The Commission suggested that upstream grouting or, if necessary, 
prestressing of the rock would prevent undesirable effects. 

The case of Bhakra dam, in this study, supports this view. 
Extensive grouting was performed in the Foandecien of Bhakra dam, 
upstream and downstream of the structure. No crack was consistent 
with the measured behaviour. However, in many cases the presence 


of a crack may be tolerable without special treatment. 
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